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ABSTRACT 

We analyze Spitzer/IRS spectra of 110 B-, A-, F-, and G-type stars with 
optically thin infrared excess in the Scorpius-Centaurus (ScoCen) OB association. 

The age of these stars ranges from 11-17 Myr. We £t the infrared excesses 
observed in these sources by Spitzer/IRS and Spitzer/MIPS to simple dust models 
according to Mie theory. We hnd that nearly all the objects in our study can be fit 
by one or two belts of dust. Dust around lower mass stars appears to be closer in 
than around higher mass stars, particularly for the warm dust component in the 
two-belt systems, suggesting mass-dependent evolution of debris disks around 
young stars. For those objects with stellar companions, all dust distances are 
consistent with trunction of the debris disk by the binary companion. The gaps 
between several of the two-belt systems can place limits on the planets that might 
lie between the belts, potentially constraining the mass and locations of planets 
that may be forming around these stars. 
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Introduction 


High contrast imaging surveys using adaptive optics (AO) enabled large telescopes are 
beginning to discover Jovian planets in nascent planetary systems. These surveys typi¬ 
cally target nearby (<100 pc), young (<300 Myr) stars because atmospheric modeling of 
Jovian planets predicts that their self-emission is bright when they are young and fades 
with time, because of th eir captured heat of formation and ongoing gravitational contrac¬ 
tion (iMarlev et ahN2007n . Debris disks are dusty disks around main sequence stars that 
are typically discovered via excess thermal infrared emission above the stellar photosphere. 
Spitzer MIPS surveys of young clusters and fi eld stars indica t e that young stare gene rally 
possess larger infrared excess than old stars flSu et al.l 120061: ICarpenter et ahl l2009all . In 
fact, t he first images of J ovian exoplanet s were made in the du sty debris disk systems, HR 
8799 flMarois et al.ll2010l) and /3 Pictoris fiLagrange et al.ll2010[) . Mid- to far-infrared excess 
was origi nally discovered toward these targets mo re than two decades ago using the IRAS 
satellite flAnmannIll984J: iRackman &: ParescdIl993l) . 


Studies of the spatial distribution of dust in debris disks can provide constraints on 
the presence of planets. For example, SED modeling of the thermal emission from the 
dust around HR 8799 indicates the presence of two dust populations gen erated by two 


2009 


separate parent b ody belts, one interior to 10 AU and another beyond 100 AU (iSu et ah 
Chen et al.ll200911 . The system of four ~10 Mj^p planets discovered using the K eck Telescope 


lies b etween these two dust populations at distances of 15, 24, 38, and 68 AU flMarois et ah 
2010l) . Similarly, SED modeling of I RAS excesses toward /3 P ictoris indicates that the 

and the Jovian planet 13 
2010). SMA observations 


dust is located at distances >20 AU (Hackman &: Paresce Il9931 


Pic b lies inside the central clearing at 10 AU fiLagrange et al 


lat has been 


of (3 Pic have further revealed a ring of millimeter-sized grains at ~94 AU t 
ident ified as the location of the main reservoir of dust-producing planetesimals (jWilner et al. 


201l|). The architecture of the {3 Pic system is consistent with creation of dust by collisions 
among parent bodies in the main belt, the larger of which spiral inward under Poynting- 
Robertson drag until they encounter (3 Pic b and are scattered out of the system. A study 
of strong silicate emis sion at 8-13 /rm of (3 Pic finds evidence for belts at 6, 16, and 30 AU 
flOkamoto et al.ll2004i) . These belts may also be sculpted by (3 Pic b, and possibly additional 
as-yet-undiscovere d planets. The infrared spectrum of the dust in (3 Pic support this picture 
flChen et al.ll2007|) . Stellar activity and starspots in young stars make the detection of planets 
by radial velocities and transits infeasible, especially in later type stars, meaning that direct 
imaging and gap characterization may be some of the best methods of finding planets in 
young stars, particularly in pre-main sequence stars. 


Simple black body modeling of ~500 debris disks observed with the Spitzer IRS and 
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MIPS at 70 /im indicates that the excess from one-third of the targets can be described 
using a single temperature black body with a median grain temperature, Tgr ~ 180 K, and 
the excess from two-thirds of the targets can be described using a two temperature blac k 
body model with median grain temperatures, Tg^ ~ 80 K and 340 K fjChen et al.N2014f) . 
In these systems, the presence of Jovian planets could naturally explain how planetesimals 
populations are (1) dynamically excited leading to collisions between parent bodies and (2) 
sculpted into rings. However, coagulation N-body simulations of ’self-stirred’ disks suggest 
that signihcantly smaller P luto-sized objects may also induce collisions between parent bodies 
flKenvon fc Bromlevll20041) . A detailed census of Jovian-mass planets in debris disks is needed 
to determine the role that Jovian planets play in exciting and sculpting parent body belts. 
Gemini South and VLT have recently commissioned GPI and SPHERE, second-generation 
coronagraphs, that are expected to take of census of planets with masses >1 Mjup around 
nearby, young stars flMacintosh et ahl 120141) . 


Prime targets for these searches will be young stars in the Scorpius-Gentaurus OB associ¬ 
ation (ScoGen). ScoGen is the closest OB association to the Sun with typical stellar distances 
of ~100 - 200 pc and contains three subgroups: Upper Scorpius (US), Upper Gentaurus Lu¬ 
pus (UGL), an d Lower Gentaurus Crux (LGG), with est imated ages of ~11 Myr, ~15 Myr, 
and ~17 Myr, fjPecaut et aLll2012l : iMamaiek et aLll2002l) respectively. Several hundred can¬ 
didate members have been identihed to date, although the association probably contains 
thousands of low-mass members. Member stars with spectral-type F and earlier have been 
i dentihed using moving group analysis of Hipparcos positions, parallaxes, and proper motions 
fide Zeeuw et al.lll999l) . while later-type members have been identih ed using youth indica¬ 
tors (i.e., high coronal X -ray activity and large lithium abundance; fjPreibisch fc Mamaiek 
2008; Slesnick et ah 20061 )1. Jovian mass planets have already been discovered in two ScoGen 
debris disk systems thus far. VLT NaCO diherential imaging at L’-band has revealed the 
presence of a 5.2 Mj „p planet at 56 AU (0.62") from HD 95086, an A8V member of LGG 
(IRameau et ahl 120131 1. Magellan AO -|- Cho2 diherential imaging at J-, Ks- and L’-bands 
has reve aled the presence o f an 11 Mjup planet at 650 AU from HD 106906, a F5 member 
of LGG ( Bailey et ah 2014 '). 


We report here the results of a study modeling the Spitzer IRS and MIPS 70 /rm 
SEDs of all of the debris disks around B- through G-type ScoGen members observed during 
the Spitzer cryogenic mission. Our scientihc goal is to better constrain the location of 
debris dust and infer the presence of planets and their orbital properties where possible. 
We list the targets for the sample, along with their spectral types, distances, and subgroup 
memberships in Table 1. For debris disks with two belts, the width of a gap between the 
belts can provide importaiit con st raints on the mass of a p lanet orbiting within the disk 
f Quillen 2006 : Chiang et ah 2009). Nesvold &: Kuchner ( 2014 ') used the 3D collisional debris 
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disk model SMACK to derive a relationship between gap width, planet mass, stellar age, 
and disk optical depth. We nse this relationship to determine which of the two-belt ScoCen 
targets are consistent with a planet on a circular orbit, and which require multiple planets 
or eccentric planet orbits. For the systems consistent with a single non-eccentric planet, we 
place an upper limit on the mass of the putative perturbing planet. 


Observations 


The infrared properties of nearby, young stars in ScoCen were systematically explored 
during the Spitzer cryogenic mission. Observers used the MIPS mid-IR photometric camera 
to search for infrared excess at 24 and/or 70 pm around ScoCen memb ers selected based on 

Hiypa r cos astrometry (de Zeeuw et al.il999 ), col or-magnitude diagrams Preibisch fc Zinneckei _ 

1 1999 ): Preibisch et ah ( 2002h and x-ra v surveys ( Walter et al.lll994 : Martfnlll998 : IPreibisch et al, 


1998 


Kunkellll999l: iKohler et al.ll2000[) . They followed-up excess targets using the IRS mid- 


IR spectrometer at 5-35 pm to search for solid-state emission features and characterize the 
shapes of the SEDs. Taken together, the Spitzer MIPS photometry indicated that approx- 
i mately one-qu a rter of the ^600 ScoCe n stars observed u s ing M IPS possess infrared excess 


1 Su et al.l 120061: ICarpenter et al.l l2009bl: IChen et al.ll201ll. 120121) . For stars with ages 10 - 
20 Myr, stellar evolutionary models suggest that stars with spectral type earlier than F are 
main sequence while stars later than F are not. The infrared excess properties of late-type 
stars in Upper Sco are consistent with gas-rich, optically thick T Tauri stars while those of 
early-type stars in Upper Sco and early- and solar-type stars in UCL and LCC are consistent 
with gas-depleted, optically thin debris disks. Because our models are only valid for optically 
thin disks, we focus only on the debris disks. 

The Spitzer IRS spectra for ScoCen members possess a diverse array of properties. For 
example, the spectrum of the F3/F5 LCC member HIP 63975 (HD 113766) shows prominent 
10 and 20 /im silicate emission features consistent with the presence of forsterite, enstatite, 
olivine, and pyroxene-rich dust, generated by the destruction of a > 300 km radius asteroid. 
Detailed modeling of the IRS spectrum suggests that the dust in this system is located in 
two cold belts located at 4 - 9 AU and 30 - 80 AU from the star, plus a warm belt at 1.8 AU 


(iLisse et al.ll2008[) . By contrast, the IRS spectra of debris disks around B- and A-type stars in 
Upper Sco reveal rising continuua without strong solid-state emission f eatures. The dust in 
these systems has been modeled using a single temperature black body flDahm fc Carpenter 
20091) . We have collected all of the IRS spectra of ScoCen US, UCL, and LCC debris disks and 
analyzed their spectra self-consistently. In Figure 1, we plot the iCs-[24] color as a function 
of J — 77 color (as a proxy for spectral type) for all of the sources observed using MIPS and 
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Fig. 1.— Color-color plot of all Sc oCen targets surveyed with Spitzer. We plot in black all 
of th e obje cts surveyed with MIPS f Su et al. 2006; Carpenter et al. 2008 . 2009b; Chen et al. 
20 111. I 2 OI 2 I) . and in red the objects whose IRS spectra are studied here. 


overlay the targets analyzed here in red. We note that the IRS spectra for 26 disk-bearing 
members of the ~11 Myr old Upper Sco have been modeled in detail by lDahm fc Carpenter 


(120091 ): for self-consistency, we independently model the spectra of all of the debris disks in 
their sample but do not reanalyze those of the primordial disks. We further note that some 
MIPS UCL and LCC excess sources were not observed using the IRS. 

We drew the calibrat ed IRS spectra for our targets from the Spitzer IRS Debris Disk 
Catalog (jChen et al.ll2014l ). Calibrated IRS low-resolution spectra typically possess absolute 
calibration uncertainties of 5% while calibrated MIP S 24 /xm fluxes typically possess cali¬ 
bration uncertainties of 2% flEngelbracht et al.l 120071) . Therefore, the spectra in the Debris 
Disk catalog are pinned to the MIPS 24 pm fluxes as reported in the Spitzer Enhanced 
Imaging Products (SEIP) Catalog to improve the absolute calibration of the data. The 
IRS Debris Disk catalog contains not only the spectra for debris disks but also their re- 


^http://irsa.ipac.caltech.edu/data/SPITZER/Enhanced/Imaging/overview.htinl 
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peatability uncertainties, ctirs.a, estimated using the difference between spectra obtained 
at two separate nod positions. Since the repeatability uncertainty can vary substantially 
from pixel-to-pixel, inconsistent with our understanding of the instrument, we averaged the 
repeatability uncertainty in quadrature over the nearest 5 points with boxcar weighting to 
smooth out anomalously low or high values. To wit, if cxo,* is the IRS repeatability error at 
Aj, then 


a 


2 

IRS,i 


i+2 

E 

j=t-2 




0,i 


( 1 ) 


One source, HIP 77911, was observed only at high resolution with IRS. For this source, 
we used the CASSIS ('Co rnell AtlaS of Spitzer/IRS Sources) optimal reduction of the data 
(ILebonteiller et al.N2nilh . The spectrum appears different from that of the other sources 
because the spectral resolution is higher and the wavelength coverage does not extend short- 
ward of 10 microns. 

We photosphere-subtracted our spe ctra using our ow' t i stell ar photosphere models. For 
stars whose MIPS data were analyzed bv IChen et al.l (1201 11.1201211 . we adopted stellar spectral 
types, effective temperatures (Teff), visual extinctions (Rv), and luminosities (L*) published 
therein. Stellar properties for sources not analyzed in these papers were t aken from the 
references indicated in Tabled! Then, we selected Kurucz model atmosphere^ (Kurucz 1979 1 
consistent with the listed effective temperatures, assuming solar abundances and surface 
gravities, log^f = 4.0. Next, we reddened the stellar model SEDs assuming the Cardelli, 
Clayton, & Mathis interstellar extinction law and Ay = 3.1E(B — V). Finally, we normalize d 
these stellar atmospheres to the MIPS 24 /rm predictions given bv IChen et ahl (120 111. l2012l) . 
with the same 3% photosphere calibration error used in those works. The total uncertainty 
in the photosphere-subtracted spectrum is 


^excess, A — 


‘^IRS,A + '^phot.A 


( 2 ) 


where cxphot = 3% x Fphot- In Figures [2all2el we show the reduced and photosphere-subtracted 
spectra of our sources. 


In Table dl we summarize the stellar properties of the 119 Sco Cen members whose 
IRS spectra ar e discussed here. The distances to the sources are taken from Hipparcos 
measurements (Ivan LeeuwenI 1200711 . Of these sources, 5 are Be stars, and are not analyzed: 


HIP 63005, HIP 67472, HIP 69618, HIP 77859, and HIP 78207. Another 4 sources are 
optically thick protoplanetary disks: HIP 56354 (HD 100453), HIP 56379 (HD 100546), HIP 


^http://www.stsci.edu/hst/observatory/cdbs/k93models.html 
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Fig. 2a.— Photosphere-subtracted spectra. The IRS spectrum is shown in gray and the 
photosphere fit is shown in cyan. The parameters for the photosphere fit for each star are 
listed in Table [U The photosphere-subtracted spectrum is plotted in black, with gray error 
bars. MIPS photometric points are indicated by magenta points, with arrows indicating 
upper limits. 
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Fig. 2b.— Continuation Figure [2al spectra of objects 
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Fig. 2c.— Continuation Figure l2al spectra of objects 
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Fig. 2d.— Continuation Figure [2al spectra of objects 
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Fig. 2e.— Continuation Figure l2al spectra of objects 
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Fig. 3.— Lir/L^ versus stellar mass. Asterisks are Lower Centaurus Crux, squares are 
Upper Centaurus Lupus, and filled triangles are Upper Scorpius. 


77157 (HT Lupi), and HIP 82747 (AK Sco. See e.g. iManoi et al.ll2006l: ISturm et al.ll2013[) . 
These sources cannot be adequately modeled using simple grain models, so we exclude them 
from our study. In Table |2] we tabulate the calculated Lir/L* based on the calibrated and 
photosphere-subtracted spectra of the remaining 110 sources. In Figure [3l we plot L/r/L* 
versus stellar mass. The downward trend versus stellar mass is contrary to what would be 
expected if disk temperatures or masses simply scale with stellar mass. This can be explained 
either by decreasing disk mass or decreasing dust temperature with increasing stellar mass, 
scenarios that will be discussed in our Results. 

In general, we found that the photosphere models were consistent with the IRS observa¬ 
tions for all of the stars in our study with the exception of HIP 56673 (HD 101088) and HIP 
78977 (HD 144548). For these two objects, this normalization of the IRS spectra produced a 
Rayleigh-Jeans power-law excess at 5-30 pm excess that could indicate the presence of a hot 
dust component with Tgr ^ 500 K. Visual spectra of HIP 56673B (HD 101088 B) show time- 


varia ble Ha emission, consistent with accretion observed toward T Tauri stars (IBitner et ah 


2010l ). Alternatively, this mismatch in the MIPS 24 pm flux and the IRS spectrum could 
indicate that these sources possess time variable excesses similar to that observed toward ID 
8 caused by stochastic grinding events (IMeng et al.ll2014l) . 


Several sources show little to no excess in the IRS data. We omit these non-excess 
sources based on two criteria: the excess significance and the normalized flux ratio. These 
quantities compare the observed versus predicted photosphere-only emission over a selected 
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wavelength regime. That is, if Th is the frequency-dependent flux with uncertainty a^, then 

CF.du 


F(z/i, z/ 2 ) = 


1^2 - 


( 3 ) 


and the weighted uncertainty is 


al du 


z/2 - Z/1 


( 4 ) 


We use the subscripts ’obs’ and ’pred’ to refer to the observed and predicted flux, respectively. 
We consider three different passbands: 8.5 — 13 pm, 21 — 26 pm, and 30 — 34 pm. The flux 
integrated over each of these passbands are F(10 pm), F(24 pm), and F(32 pm), respectively. 


The excess significance is defined to be 


X ~ (-^obs -^pred)/(c’'obs 4" ‘^pred) ^ ! 


( 5 ) 


where cxobs includes the repeatability error and a 5% normalization uncertainty added in 
quadrature, while Upred consists of a 3% normalization uncertainty. We calculate the excess 
significance over our three passbands, XiO; X 24 , and X 32 , and list their values in Table |2l 
We also list values of Xtot, which is the excess significance calculated over the entire IRS 
spectrum. HIP 56673 and HIP 78977 are listed twice : in the firs t listin g, the IRS spectra 
are normalized to the MIPS 24 micron as described in IChen et al.l (120141) . and in the second 
listing, marked by an asterisk, the spectra are normalized to the photosphere model at 5-6 
microns. When these sources are normalized to the photosphere, HIP 56673 exhibits a small 
excess, but HIP 78977 has none. 


The normalized flux ratio is adapted from ICarpenter et ahl f)2009bl) . and is defined to be 


R 


32/10 


Fobs(32pm)/Fpred(32 pm) 
Fobs (10 pm)/Fpred(10 pm)' 


( 6 ) 


A similar expression is used to calculate F 24 /io- In ICarpenter et ahl (1200911) . photometric 
fluxes were used. Here, we integrate the spectrum over the given passband, assuming 100% 
efficiency. Since R can be considered to be a ratio of the slope of the observed spectrum 
compared to the slope of the predicted spectrum, any error in the overall normalization of 
either spectrum cancels out. Therefore, the error on R is propagated from the repeatability 
error of the observed spectrum alone. 


The excess significance (y) measures the signal-to-noise of the infrared excess at each 
band pass, while R measures the shape of the excess. In some sources, even though the excess 
significance is formally low, the shape of the spectrum rises at long wavelengths, indicating 
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that there is a notable cold excess. If i? > 1, then the shape of the spectrum indicates a cold 
excess component, while R = 1 indicates a shape consistent with photospheric emission. 
To determine which spectra to exclude as non-excess sources, we use both the x and R 
measures. Sources which have Xsi X 24 , Xso; and Xtot all less than 3 and whose values of R 24/8 
and i? 3 o /8 are both less than or equal to 1 to within one a are labeled non-excess sources and 
are excluded from further analysis. A total of 13 of our sources are non-excess sources, and 
are labeled as such in Table [2l In addition, HIP 78977, when normalized to the photosphere 
model, can be considered a non-excess source. This leaves a total of 97 debris disk spectra 
that we analyze for their dust properties. 


Table 1. Stellar Properties 


HIP 

Name 

Spectral 

Distance 

TeS 

Mass 

Luminosity Ay 

Program 

Notes 



Type 

(pc) 

(K) 

{Mq) 

{Lq) 




Lower Centaurus Crux 


53524 

HD 95086 

A8HI (6) 

90.4 

7499 

1.6 

7.11 

0.000 

IRS.DISKS/2 

(18) 

55188 

HD 98363 

A2V (6) 

123.6 

8770 

1.9 

11.4 

0.244 

CCHEN2/40235 

(18) 

56354 

HD 100453 

A9Ve (14) 

121.5 

7447 

1.6 

10.4 

0.214 

IRS.DISKS/2 

(18) protoplanetary 

56379 

HD 100546 

B9Vne (6) 

96.9 

10520 

2.4 

26.6 

0.194 

IRS.DISKS/2 

(18) protoplanetary 

56673 

HD 101088 

F5IV (6) 

93.8 

6440 

2.2 

17.7 

0.142 

DEBRISH/40651 

(17) excess 

57524 

HD 102458 

F9IV (13) 

91.7 

6115 

1.2 

1.89 

0.178 

DEBRISH/40651 

(17) 

57950 

HD 103234 

F2IV/V (6) 

98.1 

6890 

1.5 

3.90 

0.066 

DEBRISH/40651 

(17) 

58220 

HD 103703 

F3V (6) 

98.9 

6740 

1.5 

3.35 

0.096 

YOUNG A/84 

(17) 

58528 

HD 104231 

F5V (6) 

110.5 

6440 

1.4 

3.73 

0.016 

DEBRISH/40651 

(17) 

58720 

HD 104600 

B9V (6) 

105.7 

11614 

2.7 

68.7 

0.024 

CCHEN2/40235 

(18) 

59282 

HD 105613 

A3V (6) 

104.2 

8551 

1.8 

11.6 

0.192 

CCHEN2/40235 

(18) 

59397 

HD 105857 

A2V (6) 

113.0 

8770 

1.9 

14.0 

0.132 

CCHEN2/40235 

(18) 

59481 

HD 105994 

F3V (7) 

113.1 

6740 

1.5 

4.09 

0.023 

WARMDISK2/50538 

(17) 

59502 

HD 106036 

A2V (6) 

100.7 

8770 

1.9 

13.8 

0.015 

CCHEN2/40235 

(18) 

59693 

HD 106389 

F6IV (7) 

137.0 

6360 

1.3 

2.32 

0.184 

DEBRISH/40651 

(17) 

59898 

HD 106797 

AOV (6) 

96.0 

9550 

2.1 

32.9 

0.022 

CCHEN2/40235 

(18) 

59960 

HD 106906 

F5V (6) 

92.1 

6440 

1.5 

5.06 

0.000 

IRS.DISKS/2 

(17) 

60183 

HD 107301 

B9V (6) 

93.9 

10814 

2.4 

36.8 

0.067 

CCHEN2/40235 

(18) 

60348 

HD 107649 

F5V (7) 

93.7 

6440 

1.4 

2.13 

0.022 

YOUNG A/84 

(17) 

60561 

HD 107947 

AOV (6) 

91.1 

9550 

2.1 

17.5 

0.000 

CCHEN2/40235 

(18) 

60710 

HD 108257 

B3Vn (2) 

137.4 

17298 

5.4 

809. 

0.070 

IRS.DISKS/2 

(18) 

61049 

HD 108857 

F7V (6) 

97.0 

6280 

1.4 

3.13 

0.143 

CCHEN2/40235 

(17) 

61087 

HD 108904 

F6V (6) 

97.5 

6360 

1.5 

4.97 

0.050 

CCHEN2/40235 

(17) 

61684 

HD 109832 

A9V (6) 

111.9 

7447 

1.6 

7.14 

0.261 

CCHEN2/40235 

(18) 

61782 

HD 110058 

AOV (7) 

107.4 

9550 

2.1 

10.2 

0.436 

IRS.DISKS/2 

(18) 



Table 1—Continued 


HIP 

Name 

Spectral 

Distance 

Teff 

Mass 

Luminosity 

Ay 

Program 

Notes 



Type 

(pc) 

(K) 

{Mq) 

{Lq) 




62134 

HD 110634 

F2V (7) 

115.6 

6890 

1.5 

3.74 

0.026 

YOUNGA/84 

(17) 

62427 

HD 111103 

F8 (1) 

142.7 

6200 

1.4 

3.17 

0.000 

DEBRISH/40651 

(17) 

62445 

HD 111170 

G4.5IVe (13) 

130.5 

5728 

1.6 

4.17 

0.660 

RUBBLE/148 

(13) 

62657 

HD 111520 

F5/6V (7) 

108.6 

6400 

1.3 

2.60 

0.028 

CCHEN/241 

(17) 

63005 

HD 112091 

B5V(e) (2) 

124.8 

16634 

5.0 

541. 

0.200 

CCHEN2/40235 

(18) classical Be 

63236 

HD 112383 

A2IV/V (6) 

110.7 

8770 

1.9 

20.0 

0.000 

CCHEN2/40235 

(18) 

63439 

HD 112810 

F3/5IV/V (7) 

143.3 

6590 

1.4 

3.52 

0.000 

DEBRISH/40651 

(17) 

63836 

HD 113524 

F6/8 (7) 

107.4 

6280 

1.3 

2.31 

0.000 

DEBRISH/40651 

(17) 

63839 

HD 113457 

AOV (6) 

99.4 

9550 

2.1 

20.3 

0.000 

CCHEN2/40235 

(18) 

63886 

HD 113556 

F2V (6) 

106.7 

6890 

1.5 

4.91 

0.024 

IRS_DISKS/2 

(17) 

63975 

HD 113766 

F3/5V (7) 

122.5 

6590 

1.9 

11.9 

0.000 

IRS_DISKS/2 

(17) 

64053 

HD 113902 

B8/9V (7) 

100.1 

11695 

2.7 

77.6 

0.070 

CCHEN2/40235 

(18) 

64184 

HD 114082 

F3V (6) 

85.5 

6740 

1.5 

3.18 

0.078 

IRS_DISKS/2 

(17) 

64877 

HD 115361 

F5V (6) 

125.0 

6440 

1.5 

5.02 

0.000 

CCHEN2/40235 

(17) 

64995 

HD 115600 

F2IV/V (6) 

110.5 

6890 

1.5 

4.79 

0.000 

IRS_DISKS/2 

(17) 

65089 

HD 115820 

A7/8V (7) 

96.5 

7656 

1.7 

4.83 

0.026 

CCHEN2/40235 

(18) 

65875 

HD 117214 

F6V (6) 

110.3 

6360 

1.6 

5.64 

0.000 

IRS_DISKS/2 

(17) 

65965 

HD 117484 

B9.5V (7) 

147.3 

10593 

2.4 

25.9 

0.100 

CCHEN/40235 

(18) 

66001 

HD 117524 

G2.5IV (13) 

152.4 

5834 

1.2 

2.14 

0.110 

RUBBLE/148 

(13) 

66068 

HD 117665 

A1/2V (7) 

147.9 

8974 

1.9 

24.4 

0.000 

GCHEN/40235 

(18) 

66566 

HD 118588 

AIV (7) 

126.4 

9204 

2.0 

14.9 

0.097 

GCHEN/40235 

(18) 

67068 

HD 119511 

F3V (7) 

91.6 

6740 

1.5 

2.70 

0.000 

WARMDISK2/50538 

(17) 

67230 

HD 119718 

F5V (6) 

131.8 

6440 

1.8 

8.67 

0.037 

CCHEN2/40235 

(17) 

Upper Centaurus Lupus 

66447 

HD 118379 

A3IV/V (7) 

121.7 

8551 

1.8 

13.3 

0.184 

TD_GTO/50485 

(18) 

67472 

HD 120324 

B2V:e (8) 

155.0 

20512 

7.3 

6.53e+03 

0.287 

CCHEN2/40235 

(18) classical Be 

67497 

HD 120326 

FOV (7) 

107.4 

7200 

1.6 

4.45 

0.158 

DEBRISH/40651 

(17) 






Table 1—Continued 


HIP 

Name 

Spectral 

Type 

Distance 

(pc) 

TeS 

(K) 

Mass 

{Mq) 

Luminosity 

{Lq) 


Program 

Notes 

67970 

HD 121189 

F3V (7) 

118.8 

6740 

1.5 

3.85 

0.070 

GGHEN2/40235 

(17) 

68080 

HD 121336 

AlVn (11) 

139.9 

9204 

2.0 

64.9 

0.082 

GGHEN2/40235 

(18) 

68781 

HD 122705 

A4V (5) 

112.9 

8279 

1.8 

8.89 

0.000 

TD.GTO/50485 

(18) 

69291 

HD 123889 

F2V (9) 

132.3 

6890 

1.5 

5.02 

0.043 

WARMDISK2/50638 

(17) 

69618 

HD 124367 

B4Vne (2) 

147.7 

16982 

5.2 

l.Ole+03 

0.409 

GGHEN2/40235 

(18) classical Be 

69720 

HD 124619 

FOV (6) 

133.3 

7200 

1.6 

5.00 

0.228 

DEBRISH/40651 

(17) 

70149 

HD 125541 

A9V (7) 

113.3 

7447 

1.6 

3.18 

0.146 

TD.GTO/50485 

(18) 

70441 

HD 126062 

AIV (7) 

110.4 

9204 

2.0 

11.4 

0.018 

GGHEN2/40235 

(18) 

70455 

HD 126135 

B8V (7) 

165.0 

11967 

2.8 

70.3 

0.096 

GOWERNER2005/20132 

(18) 

71271 

HD 127750 

AOV (7) 

175.7 

9550 

2.1 

26.9 

0.040 

TD.GTO/50485 

(18) 

71453 

HD 128207 

B8V (9) 

147.5 

13490 

3.4 

210 . 

0.010 

GOWERNER2005/20132 

(18) 

72033 

HD 129490 

F7IV/V (7) 

155.8 

6280 

1.5 

5.46 

0.297 

DEBRISH/40651 

(17) 

72070 

HD 129590 

GIV (13) 

132.6 

5945 

1.3 

2.84 

0.047 

GGHEN2/42035 

(17) 

73145 

HD 131835 

A2IV (9) 

122.7 

8770 

1.9 

10.5 

0.187 

GGHEN2/40235 

(18) 

73341 

HD 132238 

B8V (9) 

162.6 

12359 

3.0 

112 . 

0.040 

GOWERNER2005/20132 

(18) 

73666 

HD 133075 

F3IV (9) 

151.5 

6740 

2.1 

16.5 

0.364 

DEBRISH/40651 

(17) 

73990 

HD 133803 

A9V (9) 

124.8 

7447 

1.6 

8.04 

0.211 

TD.GTO/50485 

(18) 

74499 

HD 134888 

F3/5V (9) 

89.9 

6590 

1.5 

2.06 

0.058 

DEBRISH/40651 

(17) 

74752 

HD 135454 

B9.5V (7) 

173.3 

10351 

2.3 

66.2 

0.023 

GOWERNER2005/20132 

(18) 

74959 

HD 135953 

F5V (9) 

133.2 

6440 

1.3 

2.68 

0.080 

TD.GTO/50485 

(17) 

75077 

HD 136246 

AIV (9) 

131.6 

9204 

2.0 

22.7 

0.114 

GOWERNER2005/20132 

(18) 

75151 

HD 136347 

B9IVSi(SrCr) (11) 

143.3 

11641 

2.7 

68.2 

0.031 

GOWERNER2005/20132 

(18) 

75210 

HD 136482 

B8/9V (9) 

136.2 

11324 

2.6 

54.2 

0.033 

GOWERNER2005/20132 

(18) 

75304 

HD 136664 

B4V (2) 

159.2 

16711 

5.0 

1.27e+03 

0.052 

IREXT/20294 

(18) 

75491 

HD 137057 

F3V (9) 

168.6 

6740 

1.9 

9.59 

0.051 

GGHEN2/40235 

(17) 

75509 

HD 137119 

A2V (9) 

107.2 

8770 

1.9 

9.02 

0.067 

GCHEN2/40235 

(18) 

76084 

HD 138296 

F2V (9) 

142.7 

6890 

1.7 

6.64 

0.191 

DEBRISH/40651 

(17) 

76395 

HD 138923 

B8V (3) 

106.5 

11967 

2.8 

53.6 

0.041 

GOWERNER2005/20132 

(18) 

77081 

HD 140374 

G7.5IV (13) 

200.8 

5521 

1.4 

2.24 

0.150 

RUBBLE/148 

(13) 
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HIP 

Name 

Spectral 

Distance 

Teff 

Mass 

Luminosity 

Ay 

Program 

Notes 



Type 

(pc) 

(K) 

{Mq) 

(Lq) 




77157 

HT Lupi 

K3Ve (15) 

141.2 

4730 

1.1 

5.09 

1.138 

CCHEN2/40235 

(17) protoplanetary 

77315 

HD 140817 

AOV (9) 

147.3 

9550 

2.1 

41.5 

0.105 

CCHEN2/40235 

(18) 

77317 

HD 140840 

B9/A0V (9) 

125.8 

10593 

2.4 

21.5 

0.029 

CCHEN2/40235 

(18) 

77432 

HD 141011 

F5V (7) 

96.3 

6440 

1.4 

1.90 

0.000 

DEBRlSH/40651 

(17) 

77520 

HD 141254 

F3V (9) 

100.8 

6740 

1.5 

1.87 

0.194 

WARMD1SK2/50538 

(17) 

77523 

HD 141327 

B9V (9) 

195.3 

10304 

2.3 

48.3 

0.139 

TD.GTO/50485 

(18) 

77656 

HD 141521 

G51V (13) 

140.1 

5702 

1.4 

2.57 

0.520 

RUBBLE/148 

(13) 

78043 

HD 142446 

F3V (9) 

144.3 

6740 

1.5 

4.66 

0.124 

TD.GTO/50485 

(17) 

78555 

HD 143538 

FOV (9) 

106.3 

7200 

1.6 

3.55 

0.178 

GGHEN/241 

(17) 

78641 

HD 143675 

A51V/V (9) 

113.4 

8072 

1.7 

6.32 

0.049 

GGHEN2/40235 

(18) 

78756 

HD 143939 

B9H1 (10) 

144.7 

10740 

2.4 

40.6 

0.000 

TD.GTO/50485 

(18) 

79400 

HD 145357 

A5V (7) 

146.8 

8072 

1.7 

12.8 

0.342 

TD.GTO/50485 

(18) 

79516 

HD 145560 

F5V (7) 

133.7 

6440 

1.4 

3.84 

0.001 

GGHEN/241 

(17) 

79631 

HD 145880 

B9.5V (9) 

127.9 

10593 

2.4 

35.3 

0.355 

TD.GTO/50485 

(18) 

79710 

HD 145972 

FOV (7) 

127.4 

7200 

1.6 

5.78 

0.094 

TD.GTO/50485 

(17) 

79742 

HD 146181 

F6V (19) 

146.2 

6360 

1.4 

3.66 

0.000 

GGHEN/241 

(17) 

80142 

HD 147001 

B7V (7) 

137.2 

11912 

2.8 

78.5 

0.166 

TD.GTO/50485 

(18) 

80897 

HD 148657 

AOV (9) 

165.6 

9550 

2.1 

21.3 

0.864 

TD.GTO/50485 

(18) 

82154 

HD 151109 

B8V (4) 

222.7 

11298 

2.6 

117. 

0.074 

TD.GTO/50485 

(18) 

82747 

AK Sco 

F5V (9) 

102.8 

6440 

1.5 

4.72 

1.098 

lRSDlSKS/2 

(17) protoplanetary 

83159 

HD 153232 

F5V (9) 

146.6 

6440 

1.5 

4.17 

0.000 

GGHEN/241 

(17) 

Upper Scorpius 

76310 

HD 138813 

AOV (16) 

150.8 

9750 

2.1 

30.5 

0.155 

JMCARP/30091 

(20) 

77859 

HD 142184 

B2V (12) 

130.9 





JMCARP/30091 

classical Be 

77911 

HD 142315 

B9V (19) 

147.7 

10000 

2.5 

50.2 

0.202 

JMCARP/30091 

(20) 

78207 

HD 142983 

B8Ia/Iab (12) 

143.5 





JMCARP/30091 

classical Be 

78663 

HD 143811 

F5V (9) 

144.3 

6440 

1.5 

4.86 

0.093 

DEBRlSH/40651 

(17) 
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HIP 

Name 

Spectral 

Type 

Distance 

(pc) 

Feff 

(K) 

Mass 

(Mq) 

Luminosity 

(Lq) 

Ay 

Program 

Notes 

78977 

HD 144548 

F7V (15) 

116.7 

6280 

1.5 

5.02 

0.365 

CCHEN/241 

(17) excess 

78996 

HD 144587 

A9V (16) 

108.5 

8750 

1.4 

11.4 

0.967 

JMC ARP/30091 

(20) 

79054 

HD 144729 

FOV (12) 

138.9 

7200 

1.5 

5.47 

0.585 

CCHEN/241 

(17) 

79156 

HD 144981 

AOV (16) 

170.4 

9750 

1.9 

31.6 

0.574 

JMCARP/30091 

(20) 

79288 

HD 145263 

FOV (12) 

149.9 

7200 

1.6 

6.40 

0.399 

DEBRlSH/40651 

(17) 

79410 

HD 145554 

B9V (16) 

140.4 

10000 

1.9 

32.5 

0.577 

JMCARP/30091 

(20) 

79439 

HD 145631 

B9V (16) 

131.8 

10000 

1.8 

32.7 

0.666 

JMCARP/30091 

(20) 

79878 

HD 146606 

AOV (16) 

129.4 

10000 

2.1 

27.1 

0.016 

JMCARP/30091 

(20) 

79977 

HD 146897 

F2/3V (12) 

122.7 

6815 

1.5 

3.66 

0.341 

lRS.DlSKS/2 

(17) 

80024 

HD 147010 

B9H (16) 

163.4 

10500 

2.1 

73.7 

0.772 

JMCARP/30091 

(20) 

80088 

HD 147137 

A9V (19) 

139.1 

9000 

1.7 

12.7 

1.147 

JMCARP/30091 

(20) 

80320 

HD 147594 

G3IV (15) 

142.0 

5830 

1.4 

3.39 

0.031 

WARMD1SK2/50538 

(17) 

82218 

HD 151376 

F2/3V (12) 

135.7 

6815 

1.5 

4.48 

0.300 

CCHEN/241 

(17) 


Reference s. — R efere nces: (1~) Cannon fc Pickering ( 192C). (2) Hiltner et al. ( 19^, (3) Hubei ( 197^ , (4) 


Schildetal 


(197lli 


iGlaspevI ( 1972h. f6 i iHouk fc CowlevI (Il975ll. I7i iHouk (Il978h . (8) Morgan et al.l (11 978^. ('9 ') iHouk (1198211 . ('101 iGahm et al 


( 198,‘f[l. fll) Corballvl ( 1984li. fl2) Honk fc Smith-Moore ( 1988), (13 ') iMamaiek et al.l (2002 ) (14) Vieira et al. ( 20031). 115 ) 


Torres et al 


1201 IChen et al.l (120141) 


( 20061). d Ol TPreibisch fc MamaiekI (I2OO8I) . (^17) iHen et all d^ll) . llSl IChen et all (l2012(l . (19) 


Pecaut et al 


(I2OI2I) . 









































































































Table 2. Inferred disk properties 


HIP 

HD name Lir/L* 

excess significance (x) 

^24/10 

t?32/10 

Omin note 



Xio X24 X32 total 



(/im) 


Lower Centaurus Crux 


53524 

HD 95086 

l.lOe-03 

1.85 

7.10 

9.75 

2.13 

3.21 ±0.30 

12.53 ± 1.08 

1.8 

55188 

HD 98363 

9.56e-04 

6.47 

16.55 

7.98 

5.19 

9.97 ±0.31 

23.44 ±2.65 

2.3 

56673 

HD 101088 

5.50e-04 

2.86 

4.37 

2.90 

2.91 

1.12 ±0.02 

1.11 ±0.07 

3.2 

56673* 

HD 101088 

4.42e-05 

0.23 

2.00 

1.30 

0.28 

1.12 ±0.02 

1.11 ±0.07 

3.2 

57524 

HD 102458 

2.91e-04 

1.27 

2.56 

2.13 

1.19 

1.40 ±0.19 

1.79 ±0.41 

0.4 

57950 

HD 103234 

1.44e-04 

0.96 

7.04 

4.65 

0.97 

1.76 ±0.07 

2.45 ±0.30 

1.1 

58220 

HD 103703 

6.91e-04 

4.40 

6.30 

4.24 

3.50 

2.19 ±0.23 

3.17 ±0.60 

0.9 

58528 

HD 104231 

2.35e-04 

2.47 

7.48 

3.55 

1.15 

1.98 ±0.12 

2.57 ±0.47 

1.1 

58720 

HD 104600 

1.09e-04 

3.25 

15.10 

17.02 

2.44 

3.72 ±0.06 

8.04 ±0.16 

8.6 

59282 

HD 105613 

8.07e-05 

2.40 

9.26 

4.78 

0.94 

1.91 ±0.06 

2.58 ±0.33 

2.5 

59397 

HD 105857 

1.74e-04 

2.32 

12.06 

9.20 

1.95 

3.01 ±0.13 

4.39 ±0.34 

2.8 

59481 

HD 105994 

8.28e-05 

0.17 

1.86 

1.16 

0.42 

1.19 ±0.09 

1.39 ±0.34 

1.1 

59502 

HD 106036 

3.79e-04 

5.49 

13.80 

5.17 

3.38 

4.71 ±0.18 

8.45 ± 1.45 

2.8 

59693 

HD 106389 

3.83e-04 

3.48 

4.58 

2.27 

1.57 

1.67 ±0.18 

1.87 ±0.47 

0.7 

59898 

HD 106797 

2.21e-04 

3.85 

16.05 

17.48 

3.11 

4.31 ±0.06 

8.23 ±0.14 

5.5 

59960 

HD 106906 

1.27e-03 

0.44 

13.60 

16.09 

2.49 

7.11 ±0.35 

23.54 ± 1.04 

1.4 

60183 

HD 107301 

9.65e-05 

1.60 

15.08 

14.85 

1.74 

4.06 ±0.14 

8.46 ± 0.40 

5.4 

60348 

HD 107649 

2.27e-04 

1.50 

6.69 

3.92 

1.07 

1.83 ±0.11 

3.03 ±0.51 

0.4 

60561 

HD 107947 

1.13e-04 

2.33 

11.31 

10.33 

1.49 

2.44 ±0.08 

4.33 ±0.27 

3.1 

60710 

HD 108257 

-2.13e-06 

-0.43 

-0.11 

0.25 

-0.21 

1.02 ±0.07 

1.11 ±0.31 

44.7 

61049 

HD 108857 

4.36e-04 

3.46 

12.94 

10.42 

1.86 

3.13 ±0.09 

4.16 ±0.25 

0.9 

61087 

HD 108904 

5.24e-04 

2.43 

5.69 

1.95 

1.88 

2.93 ±0.35 

3.58 ± 1.38 

1.4 

61684 

HD 109832 

4.16e-04 

1.14 

13.01 

15.03 

1.79 

3.93 ±0.13 

10.36 ±0.39 

1.8 

61782 

HD 110058 

1.58e-03 

3.69 

17.41 

18.00 

6.63 

23.13 ± 1.24 

72.52 ±3.83 

1.9 

62134 

HD 110634 

4.61e-05 

0.08 

1.59 

0.27 

0.15 

1.59 ±0.36 

1.34 ± 1.24 

1.0 

62427 

HD 111103 

2.46e-04 

-0.57 

7.19 

5.17 

0.57 

2.76 ±0.22 

6.38 ± 1.02 

0.9 

62445 

HD 111170 

-1.57e-04 

-0.12 

0.31 

-0.01 

-0.58 

1.04 ±0.10 

1.01 ±0.39 

1.1 

62657 

HD 111520 

1.96e-03 

1.81 

14.82 

17.50 

3.49 

5.97 ±0.23 

20.37 ±0.70 

0.8 

63236 

HD 112383 

1.54e-04 

2.57 

12.60 

10.61 

1.73 

2.71 ±0.12 

3.76 ±0.26 

3.9 

63439 

HD 112810 

8.00e-04 

1.38 

6.25 

10.92 

1.18 

2.00 ±0.16 

6.48 ±0.47 

1.0 

63836 

HD 113524 

1.26e-04 

1.00 

4.34 

2.42 

0.47 

1.46 ±0.09 

2.56 ±0.66 

0.7 

63839 

HD 113457 

2.74e-04 

5.07 

15.37 

16.07 

3.34 

4.13 ±0.14 

7.12 ±0.26 

3.6 


A ^ excess 
* 


no excess 


no excess 
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HIP 

HD name 

Lir/L* 

excess significance (x) 

Xio X24 X32 total 

^24/10 

^32/10 

^min 

(^m) 

note 

63886 

HD 113556 

5.15e-04 

0.84 

4.68 

4.88 

0.75 

2.05 ±0.22 

6.74 ± 1.15 

1.3 


63975 

HD 113766 

2.23e-02 

17.40 

18.78 

14.00 

13.12 

4.00 ±0.11 

4.93 ±0.27 

2.5 


64053 

HD 113902 

4.17e-05 

2.27 

9.93 

8.92 

1.28 

1.89 ±0.04 

2.54 ±0.14 

9.6 


64184 

HD 114082 

3.01e-03 

1.83 

13.19 

9.41 

3.82 

19.43 ± 1.17 

56.92 ±5.46 

0.8 


64877 

HD 115361 

3.09e-04 

1.22 

5.85 

1.86 

1.04 

2.59 ±0.25 

5.43 ±2.41 

1.4 


64995 

HD 115600 

1.87e-03 

1.57 

12.22 

8.98 

3.49 

12.76 ±0.88 

40.36 ±4.16 

1.3 


65089 

HD 115820 

1.52e-04 

-0.32 

9.60 

6.66 

1.21 

2.46 ±0.09 

3.65 ±0.36 

1.2 


65875 

HD 117214 

2.53e-03 

2.14 

13.25 

9.56 

4.05 

13.32 ±0.72 

37.70 ±3.41 

1.5 


65965 

HD 117484 

2.37e-04 

3.10 

14.72 

14.76 

2.92 

4.93 ±0.18 

13.81 ±0.64 

3.9 


66001 

HD 117524 

7.98e-05 

1.31 

0.88 

-0.01 

0.28 

1.04 ±0.13 

0.91 ±0.42 

0.6 

no excess 

66068 

HD 117665 

1.94e-04 

3.18 

13.91 

11.39 

2.17 

3.45 ±0.10 

5.71 ±0.34 

4.6 


66566 

HD 118588 

2.24e-04 

3.00 

14.26 

11.26 

2.69 

4.32 ±0.16 

8.16 ±0.55 

2.9 


67068 

HD 119511 

5.58e-06 

0.28 

3.79 

2.30 

0.03 

1.33 ±0.07 

1.52 ±0.23 

0.7 


67230 

HD 119718 

3.39e-04 

1.51 

4.08 

3.30 

1.05 

2.99 ±0.49 

6.06 ± 1.56 

2.0 


Upper Centaurus Lupus 

66447 

HD 118379 

1.37e-04 

0.45 

10.08 

10.05 

0.80 

2.44 ±0.10 

6.29 ±0.46 

2.8 


67497 

HD 120326 

1.50e-03 

1.31 

15.53 

15.65 

3.50 

11.52 ±0.70 

31.39 ±2.00 

1.1 


67970 

HD 121189 

5.25e-04 

2.13 

11.04 

12.42 

2.10 

3.93 ±0.23 

8.62 ±0.52 

1.0 


68080 

HD 121336 

6.10e-05 

1.27 

8.88 

9.60 

0.83 

1.81 ±0.04 

2.83 ±0.14 

10.9 


68781 

HD 122705 

8.87e-05 

1.93 

7.62 

4.35 

0.96 

1.72 ±0.06 

2.18 ±0.27 

2.0 


69291 

HD 123889 

6.89e-05 

0.62 

2.22 

2.62 

0.38 

1.14 ±0.06 

1.71 ±0.28 

1.4 


69720 

HD 124619 

1.51e-04 

1.97 

5.28 

3.49 

0.98 

1.59 ±0.11 

2.18 ±0.36 

1.3 


70149 

HD 125541 

2.06e-04 

1.57 

9.20 

7.94 

1.22 

3.28 ±0.21 

6.80 ±0.68 

0.8 


70441 

HD 126062 

2.72e-04 

1.99 

9.09 

9.59 

1.57 

2.68 ±0.15 

6.52 ±0.51 

2.2 


70455 

HD 126135 

3.28e-05 

1.35 

7.62 

6.47 

1.01 

2.08 ±0.12 

3.25 ±0.33 

8.5 


71271 

HD 127750 

1.62e-04 

1.32 

10.74 

12.09 

1.10 

2.51 ±0.08 

7.01 ±0.37 

4.6 


71453 

HD 128207 

l.Ole-05 

0.34 

4.17 

5.05 

0.52 

1.29 ±0.04 

1.53 ±0.08 

19.5 


72033 

HD 129490 

-1.06e-04 

-0.20 

0.45 

0.64 

-0.48 

1.06 ±0.09 

1.23 ±0.34 

1.5 

no excess 

72070 

HD 129590 

5.22e-03 

1.15 

18.03 

18.66 

5.45 

15.21 ±0.43 

60.57 ± 1.80 

0.9 


73145 

HD 131835 

2.27e-03 

5.10 

17.74 

17.89 

6.89 

14.32 ±0.69 

49.81 ±2.48 

2.2 


73341 

HD 132238 

2.64e-05 

1.63 

7.60 

6.90 

0.93 

2.03 ±0.11 

3.06 ±0.28 

12.2 


73666 

HD 133075 

2.36e-05 

-0.13 

2.35 

2.85 

0.13 

1.25 ±0.08 

1.99 ±0.33 

3.1 


73990 

HD 133803 

3.62e-04 

2.95 

14.05 

12.14 

2.30 

4.12 ±0.16 

7.35 ± 0.44 

2.0 


74499 

HD 134888 

8.25e-04 

1.10 

8.84 

12.42 

1.27 

3.14 ±0.22 

10.34 ±0.67 

0.3 
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HIP 

HD name 


excess significance (x) 

Xio X24 X32 total 

^24/10 

R32/10 

^min 

note 

74752 

HD 135454 

7.22e-05 

1.70 

4.89 

5.94 

1.61 

1.29 ±0.08 

1.61 ±0.11 

9.7 


74959 

HD 135953 

6.31e-04 

0.80 

4.85 

10.84 

0.85 

1.60 ±0.11 

5.59 ±0.35 

0.8 


75077 

HD 136246 

5.09e-05 

0.57 

3.55 

6.70 

0.41 

1.31 ±0.07 

2.07 ±0.13 

4.1 


75151 

HD 136347 

5.01e-06 

0.49 

0.58 

0.21 

0.18 

1.01 ±0.04 

0.99 ±0.10 

8.5 

no excess 

75210 

HD 136482 

6.14e-05 

2.53 

11.56 

10.42 

1.62 

2.78 ±0.10 

5.03 ±0.32 

7.1 


75304 

HD 136664 

-3.90e-06 

-0.26 

0.32 

0.21 

-0.42 

1.04 ±0.04 

1.04 ±0.05 

73.8 

no excess 

75491 

HD 137057 

3.12e-04 

0.67 

11.41 

10.90 

1.35 

3.45 ±0.16 

7.43 ±0.51 

2.0 


75509 

HD 137119 

1.92e-04 

3.38 

9.02 

4.78 

2.13 

2.77 ±0.17 

4.03 ±0.64 

1.9 


76084 

HD 138296 

-7.25e-05 

-0.31 

0.50 

0.10 

-0.44 

1.07 ±0.09 

1.06 ±0.35 

1.6 

no excess 

76395 

HD 138923 

8.42e-05 

4.39 

12.14 

12.72 

2.77 

2.48 ±0.07 

3.61 ±0.14 

6.6 


77081 

HD 140374 

3.63e-04 

1.36 

-0.52 

-0.45 

0.96 

0.85 ±0.11 

0.68 ±0.47 

0.1 

no excess 

77315 

HD 140817 

1.54e-04 

2.69 

12.88 

13.92 

2.15 

3.39 ±0.11 

6.34 ±0.26 

6.9 


77317 

HD 140840 

1.64e-04 

1.82 

10.72 

14.39 

1.79 

3.63 ±0.19 

10.90 ±0.48 

3.3 


77432 

HD 141011 

3.19e-04 

1.65 

3.85 

3.14 

1.30 

1.74 ±0.21 

3.26 ±0.76 

0.2 


77520 

HD 141254 

9.31e-05 

1.18 

2.47 

2.58 

0.47 

1.28 ±0.13 

2.07 ±0.43 

0.2 


77523 

HD 141327 

4.59e-05 

1.86 

6.65 

4.79 

0.92 

1.77 ±0.10 

2.87 ±0.39 

7.2 


77656 

HD 141521 

-1.73e-04 

-0.41 

-0.06 

-0.03 

-0.61 

1.02 ±0.15 

1.00 ±0.33 

0.7 

no excess 

78043 

HD 142446 

7.00e-04 

1.60 

5.90 

9.73 

1.46 

2.54 ±0.25 

7.80 ±0.60 

1.3 


78555 

HD 143538 

-7.40e-05 

-0.52 

0.52 

0.34 

-0.21 

1.29 ±0.47 

1.55 ± 1.51 

0.9 

no excess 

78641 

HD 143675 

5.77e-04 

4.11 

15.67 

14.25 

3.61 

6.15 ±0.18 

14.16 ±0.66 

1.5 


78756 

HD 143939 

2.17e-05 

0.82 

5.42 

4.14 

0.52 

1.73 ±0.11 

2.53 ±0.36 

5.9 


79400 

HD 145357 

l.lle-04 

1.47 

5.35 

8.60 

0.95 

1.75 ±0.16 

2.63 ±0.20 

2.9 


79516 

HD 145560 

2.85e-03 

0.59 

14.18 

15.53 

2.82 

8.20 ±0.56 

33.35 ±2.30 

1.1 


79631 

HD 145880 

2.43e-04 

1.20 

15.19 

17.19 

1.73 

5.18 ±0.26 

17.28 ±0.88 

5.2 


79710 

HD 145972 

2.14e-04 

2.39 

7.07 

3.32 

1.26 

2.42 ±0.19 

3.27 ±0.72 

1.5 


79742 

HD 146181 

2.59e-03 

1.24 

8.71 

14.00 

2.62 

8.12 ±0.98 

32.95 ±3.14 

1.1 


80142 

HD 147001 

-1.42e-05 

-0.92 

-0.03 

0.20 

-0.47 

1.06 ±0.27 

1.19 ±0.54 

9.4 

no excess 

80897 

HD 148657 

3.71e-04 

3.05 

15.32 

16.48 

3.81 

5.96 ±0.25 

17.90 ±0.79 

3.7 


82154 

HD 151109 

8.35e-05 

1.94 

13.05 

12.35 

1.75 

3.73 ±0.12 

9.06 ±0.51 

14.6 


83159 

HD 153232 

5.23e-04 

0.32 

0.46 

0.34 

0.39 

1.81 ±1.80 

3.91 ±9.26 

1.1 

no excess 

Upper Scorpius 

76310 

HD 138813 

9.50e-04 

4.36 

18.00 

17.65 

6.55 

12.95 ±0.39 

39.23 ± 1.38 

5.2 


77911 

HD 142315 

3.52e-04 

0.11 

2.58 

1.75 

0.75 

8.83 ±4.33 

25.22 ± 16.45 

7.0 


78663 

HD 143811 

3.30e-05 

-0.15 

2.48 

2.20 

0.06 

1.34 ±0.12 

2.18 ±0.53 

1.3 
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HIP 

HD name 

Dir/D* 

excess significance (x) 

Xio X24 X32 total 

.^24/10 

-R32/10 

^min 

(^m) 

note 

78977 

HD 144548 

1.33e-03 

4.72 

1.60 

1.38 

4.52 

0.97 ±0.15 

1.50 ±0.53 

1.4 

excess 

78977* 

HD 144548 

9.55e-05 

0.21 

-0.13 

0.90 

0.10 

0.97 ±0.15 

1.50 ±0.53 

1.4 

no excess* 

78996 

HD 144587 

3.24e-04 

4.70 

13.91 

11.09 

3.37 

2.93 ±0.10 

4.80 ±0.31 

3.1 


79054 

HD 144729 

9.49e-05 

1.15 

1.09 

0.55 

0.51 

1.33 ±0.37 

1.50 ± 1.07 

1.5 

no excess 

79156 

HD 144981 

2.15e-04 

4.66 

10.36 

8.08 

3.09 

2.16 ±0.10 

3.27 ±0.28 

5.9 


79288 

HD 145263 

1.35e-02 

18.25 

19.36 

19.13 

16.29 

5.84 ±0.08 

8.87 ±0.15 

1.6 


79410 

HD 145554 

1.76e-04 

4.38 

12.65 

10.87 

2.68 

2.86 ±0.11 

4.95 ±0.32 

6.0 


79439 

HD 145631 

8.33e-05 

1.16 

6.31 

5.77 

1.20 

1.72 ±0.11 

2.47 ±0.26 

6.4 


79878 

HD 146606 

7.98e-05 

1.72 

6.34 

5.92 

1.01 

2.47 ±0.22 

4.50 ±0.59 

4.6 


79977 

HD 146897 

5.21e-03 

1.71 

15.88 

15.25 

4.94 

24.15 ± 1.87 

89.38 ±7.21 

1.0 


80024 

HD 147010 

1.17e-04 

1.70 

7.70 

7.76 

2.12 

3.54 ±0.30 

8.22 ±0.86 

11.6 


80088 

HD 147137 

3.98e-04 

2.41 

11.77 

14.09 

2.90 

3.44 ±0.16 

7.91 ±0.36 

2.9 


80320 

HD 147594 

1.37e-04 

1.09 

4.65 

3.22 

0.45 

1.44 ±0.09 

2.22 ±0.38 

1.0 


82218 

HD 151376 

2.77e-04 

0.52 

3.21 

2.78 

0.55 

2.07 ±0.35 

4.75 ± 1.39 

1.2 



*HIP 56673 and HIP 78977 are listed twice, the marked listing indicating that the spectrum was normalized to the 
photosphere model rather than the MIPS 24 micron measurement. 







24 


3. Debris Disk Modeling 

Our modeling of the dust in the debris disks under study goes beyond a simple blackbody 
model. This is because the spectral range covered by IRS includes various silicate features. 
In order to model the IRS spectra in better detail, we include grain properties such as size, 
temperature, and composition and generate model spectra using Mie theory. 


3.1. Optical Properties of Grains 


We assume that the dust is composed of amorphous silic ates of olivine and pyro xene 
composition, the o ptical constants for which are adopted from iDmschner et al.l (119951 ) and 
Jaeger et ah (1994), assuming a Mg : Fe ratio of 1 for both specie® 


In order to calculate equilibrium temperatures, we need the optical constants at short 
wave lengths as wel l. For A < 2/im, we use the optical constants for astronomical silicates 
from lDraine fc Leel (119841) . 


The minimum grain size, Umin, is estimated by assuming that radiation pressure removes 
the smallest grains if (3{= Rrad/Rgrav) > 0.5, so that 


> 


6L^(Qpr(a)) 

IQiiGM^cps 


(7) 


(lArtvmowica 119881) . where L* and M* are the stellar luminosity and mass, ps is the density 
of an individual grain, and {Qpr(a)) is the spectrum-averaged radiation pressure coupling 
coefficient, given by (Qpr(a)) = / Qpr(a, \)F\ d X/( f F\ d\) . Thes e values are tabulated in 
Chen et al.l (120111 ) for F- and G-type stars, and IChen et ahl (120121 ) for B- and A-type stars. 
In Table [1] we list the stellar mass and luminosity assumed for modeling the dust for each 
source, and in Table [21 we list the estimated minimum grain size. 


Chen et ah! (1201 ll. l2012l) estimated the color temperature of the dust from the the ratio 


of 24 to 70 micron MIPS photometric excess, and the grain distance calculated assuming 
the grains had a temperature equal to the color temperature and a size equal to an average 
grain size of (a) = 5amin/3. The IRS spectra provide more detailed information on the 
temperature, size, and composition of the dust grains than the MIPS photometry. In this 
work we allow the grain size to be a free parameter. The distribution of grain sizes depends 


^ Optical constants for all silicate species are available at http://www.astro.uni- 
j ena. de/Laboratory / O CDB / newsilicates. html. 
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on the model type implemented, as described in 1 13.31 and amin is treated as a lower limit on 
grain size. 

For compntational simplicity, we nse Mie theory to calcnlate the optical constants for 
scattering and absorption of light on particles of different sizes. The main element missing 



a), the extinction and scattering efficiencies, respectively, as a fnnction of wavelength 
A and grain radius a. Then the absorption efficiency is Scattered stellar 

light does not contribute signihcantly at the Spitzer IRS wavelengths, so we consider only 
the thermal component of emission in modeling the spectra. 


3.2. Grain Temperatures and Distances 

The dust grains are assumed to be heated by stellar irradiation and in thermal equilib¬ 
rium. The amount of radiation absorbed by a grain is the incident stellar flux modihed by 
the absorption efficiency Qabs scaled by the cross-section of the grain. The absorption effi¬ 
ciencies of olivine and pyroxene are and respectively. The fractional composition 
of olivine is fo, so the fractional composition of pyroxene is (1 — /o). The stellar flux at a 
distance r from the star is = 7 i{R^/ryiu,*(Tes). Then the total power absorbed by a 
grain is 

Pabs = ^ R,.iRs)diy (8) 

where R,* is the stellar spectrum, for which we use the best £t Kurucz models. 

The total emergent power of a grain of radius a and temperature Tgr is 

poo 

Femit= / A7r‘^a‘^Q^^^B,{T^,)diy (9) 

Jo 

where By is the Planck function. 

To hnd the equilibrium temperature of the grain, we set Remit = Rabs and hnd that 

/ oo p2 poo 

diy = ^j du (10) 


http://www-atm.physics.ox.ac.uk/code/mie/index.html 
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The absorption efficiency, Qahs, is itself a function of a and z/, and also depends on the 
composition, as determined from Mie theory and the optical constants of the constituents, 
namely oliving and pyroxene. We set equal to the the absorption efficiencies of olivine 
oi’ pyroxene depending on the grain composition. We assume a segregated 

spheres distribution for the composition of grains where each grain is either pure olivine 
or pure pyroxene, and fo and fp are the mass fraction of grains consisting of olivine or 
pyroxene, respectively. The distance of the grains derived from the grain temperatures is 
sensitive to the composition. For instance, highly reflective grains, such as ices, absorb 
energy less efficiently, therefore for the same equilibrium temperature they will be at smaller 
stellocentric distances compared to more absorptive grains. For olivine and pyroxene, we 
hnd that the differences in distances are not signihcant. We also hnd that including a wider 
range of compositions, such as crystalline silicates, does not signihcantly improve hts to the 
mid-IR spectra enough to justify the additional free parameters. Thus, the scope of this 
work is limited to olivine and pyroxene compositions only. 

To facilitate calculations of equilibrium temperatures, we tabulate values of r versus Tgr 
as a function of a and composition and interpolate on the grid. The equilibrium temperature 
of a grain of a given size is a proxy for its distance from the star. 


We assume that the disks are optically thin, so the total spectrum of the disk is the 
summation of the emission of all the grains in the disk. We dehne n{a, r) to be the number 
density of grains between distance r and r + Sr and grain size a and a + Sa grain size distri¬ 
bution as a function of a and r so that the total number of grains is N = f f n(a, r) 27ir drda 
and the total mass of the disk is 


M 


47ra^ 

3 


Pdn{a,r) 2nr drda 


( 11 ) 


where pd is the bulk density of the dust grains (3.3 g/cm^). The emitted spectrum of a single 
grain of radius a at a distance r from the star is 

Pp, Qabs(®) o, r )] 

where d is the distance between the observer and the star. Thus, the total integrated 
spectrum for a distribution of particles in a disk is 

f B^Q!,hs(ayX)B^[Tf^(T^jf,a,r)]n{a,r)27rrdrda. 


( 12 ) 
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3.3. Fitting the Spectra 

We fit our spectra to two different dust distribution models, allowing grain size, temper¬ 
ature/distance, and composition to be free parameters in order to understand the properties 
of the debris disks. The model types were 1) a single uniform grain size with a single tem¬ 
perature and composition, and 2) a two grain model where each population has a uniform 
grain size, temperature, and composition. In each case, the best htting size and temperature 
parameters are determined by minimizing the reduced xt- This was implemented using the 
IDL routine MPFITFUN. The detailed description of each model follows, including a listing 
of the free parameters for each model type. 

1. Single grain model: This model assumes a population of grains of a uniform single 
size and temperature, as if they were distributed in a ring (or a shell) of uniform 
radius around the star. The free parameters for this model are the grain radius a, 
temperature Tg^, and composition (/o,/p). The total number of particles N sets the 
overall normalization of the spectrum. Since the grain size distribution and stellocentric 
distance are delta functions, then the brightness of this model disk assuming a distance 
of d is 

Trr/^ 

A) + A)] (13) 

This mass of this disk is then 

A'jrn^ 

Mdisk = (14) 

Our assumption is that all the grains are at the same temperature. If the absorp¬ 
tion efficiencies were to differ greatly between olivine and pyroxene, then the grains 
would not be co-located at the same distance from the star. For the grain sizes and 
temperature ranges explored, the absorption efficiencies are similar enough that the 
dust grains are effectively co-located. For typical temperatures and grain sizes in our 
hts, the calculated distances between pure olivine and pure pyroxene grains differs by 
~ 10 %. 

2. Two-grain model: 

This model is a superposition of two uniform grain populations, each with its own 
independent grain size, temperature, and composition. That is, one ring of grains with 
grain radius ai, temperature Ti, composition (/oi,/pi), and number iVi, plus a second 
ring with parameters 02 , T 2 , (/ 02 , fp 2 ), and W- Then the brightness of this disk model 

is 

7rn‘^ 

F, = [UiQT^a,, u) + u)] {a,, u)B,{T,) 
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(15) 


This mass of this disk is then 



Pd 


(16) 


4. Results 


For each ht, the grain temperature, grain size, and amorphous silicate composition 
are allowed to be free parameters. That is, the single grain model has 4 free parameters 
(temperature, grain size, olivine:pyroxene ratio, and total mass), and the two-grain model 
has 8 free parameters. Once we have carried out a least-squares fit to the excess spectrum 
of the single-grain and two-grain population models, as described in section 13.31 we need to 
determine which model best describes each object. If the reduced xl value is less than 2, 
then we declare the £t to be reasonable. The resulting hts for the 96 sources analyzed are 
shown in Tables |3}l5l These tables also list the derived stellocentric distance of each grain 
population, rgr, which is calculated from the htted grain properties and assumed stellar 
properties, according to Equation (ITnl) . In Figures HaMel we show these fits to the infrared 
excess spectra. In each plot, the photosphere-subtracted spectrum is plotted, along with the 
best fit single-grain and two-grain models to the excess spectrum. 

If the single grain model is a reasonable fit for a given object, then that object is 
considered a single-belt debris system. If the single grain model is not a reasonable fit for 
the object but the two-grain model is, then we consider that object a two-belt debris system. 
Upon visual inspection, we found 10 objects whose formal xl values for a single grain model 
are less than 2, but have much better hts to the two-grain model, and are best explained 
by noisiness of the spectra. We re-categorize these objects as two-belt systems. Similarly, 
5 objects have formal xl greater than 2 for a single grain model, but are not significantly 
better ht by a two grain model, and we re-classify these as single belt systems. An additional 
11 systems have xl for both hts larger than 2, but are mostly well-ht by a two grain model, 
and we re-categorize them as two belt systems. Most of the objects in the latter two groups 
have mismatches at the short wavelength part of the spectrum, which is most ahected by 
the normalization of the stellar photosphere. HIP 55188 (Figure 14a)) is one of the systems 
for which neither model produces a good ht, although it is clear that the two-grain model 
produceds a much better ht than the single-grain model for this particular case. For each 
ht, the grain temperature, grain size, and amorphous silicate composition are allowed to be 
free parameters. (Note: HIP 55188 is well-ht if we include crystalline silicates.) In total, we 
have 48 objects that are single-belt systems and 44 that are two-belt systems. 
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Fig. 4a.— Best fit single grain and two-grain models to photosphere-subtracted infrared spectra. 
Solid lines indicate a good (x^ < 2) fit, while dotted lines indicated a poor fit. The single grain fit 
is plotted in magenta, and the two-grain fit is plotted in green, with the two components plotted 
as dashed and dot-dashed lines. The lower panel of each plot shows the residuals after the fit. 
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Fig. 4b.— Continuation Figure Hal 
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Fig. 4c.— Continuation Figure Sal 
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Fig. 4d.— Continuation Figure Hal 
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Fig. 4e.— Continuation Figure Hal 





















































































Table 3. Fits to single grain model 


HIP ID 

xl 

Tgr 

(K) 

Ogr 

(/xm) 

mass 

(Afnioon) 

fo 

rgr 

(AU) 

tHIP 56673 

0.93 

267±12 

7.59±0.83 

1.75e-05 

O.OOiO.OO 

3.69±0.33 

*HIP 56673 

2.19 

1480±40 

198±8 

9.38e-05 

O.OOiO.OO 

0.151±0.02 

HIP 57524 

1.41 

246±4 

85±4.3 

2.12e-04 

O.OOiO.OO 

1.73±0.12 

HIP 57950 

1.14 

203±2 

8.07±0.25 

4.39e-05 

O.OOiO.OO 

2.93±0.36 

HIP 58528 

1.39 

274±1 

4.68±0.08 

1.77e-05 

O.OOiO.OO 

1.63±0.17 

HIP 59282 

1.50 

274±2 

6.37±0.12 

2.94e-05 

O.OOiO.OO 

2.92±0.26 

HIP 59397 

1.27 

216±1 

5.24±0.05 

1.07e-04 

O.OOiO.OO 

5.15±0.43 

HIP 59481 

0.47 

176±4 

4.98±0.42 

1.03e-05 

O.OOiO.OO 

4.1±0.55 

HIP 59693 

1.11 

438±4 

4.28±0.1 

3.41e-06 

0.54±0.11 

0.618±0.07 

HIP 59960 

1.29 

108±0.3 

11.2±0.1 

1.20e-02 

l.OOiO.OO 

12.3±3.4 

HIP 60348 

1.95 

182±2 

8.87±0.23 

5.64e-05 

O.OOiO.OO 

2.65±0.13 

HIP 61049 

0.92 

250±0.5 

4.9±0.03 

5.47e-05 

O.OOiO.OO 

1.74±0.5 

HIP 61087 

1.85 

208±1 

3.27±0.07 

7.75e-05 

0.13±0.11 

3.49±0.38 

HIP 62134 

0.60 

175±6 

1030±160 

1.02e-02 

l.OOiO.OO 

5.3±1 

HIP 62427 

0.91 

112±1 

37.2±2 

7.02e-03 

l.OOiO.OO 

10±1.4 

HIP 63236 

1.53 

235±1 

5.22±0.03 

l.OOe-04 

O.OOiO.OO 

5.2±0.61 

HIP 63836 

1.26 

227±4 

8.92±0.36 

1.79e-05 

O.OOiO.OO 

1.81±0.3 

HIP 64053 

0.50 

294±1 

9.52±0.1 

1.08e-04 

O.OOiO.OO 

6.61±0.87 

HIP 64877 

1.00 

167±2 

8.85±0.35 

2.94e-04 

O.OOiO.OO 

4.83±0.48 

HIP 64995 

1.24 

114±1 

9.48±0.29 

1.14e-02 

0.97±0.28 

10.7±1.4 

*HIP 65089 

2.13 

93.1±0.7 

0.102±0.01 

1.70e-03 

0.59±0.07 

131±15 

HIP 66447 

0.83 

114±1 

10±0.3 

2.04e-03 

0.84±0.22 

17.7±2 

HIP 67068 

0.63 

226±8 

4.48±0.44 

4.62e-06 

0.25±0.70 

2.07±0.36 

HIP 67497 

1.53 

124±0.4 

15.5±0.2 

9.96e-03 

l.OOiO.OO 

8.82±1 

HIP 68781 

1.06 

276±2 

4.9±0.08 

1.36e-05 

O.OOiO.OO 

2.56±0.35 

HIP 69291 

1.09 

102±2 

5.89±0.43 

1.71e-04 

l.OOiO.OO 

15.3±1.7 

HIP 69720 

1.36 

285±4 

7.92±0.27 

1.63e-05 

O.OOiO.OO 

1.77±0.26 

HIP 70455 

0.94 

216±2 

8.93±0.19 

1.95e-04 

O.OOiO.OO 

11.5±1.3 

HIP 71453 

0.43 

224±6 

119±12 

2.18e-03 

0.80±0.71 

21.1±4.1 

*HIP 72070 

5.98 

93.7±0.1 

16.2±0.1 

7.98e-02 

l.OOiO.OO 

12.2±2.3 

HIP 73341 

0.44 

242±2 

9.57±0.19 

2.09e-04 

O.OOiO.OO 

11.9±1.6 

HIP 73666 

0.83 

132±4 

14.8±1.5 

7.65e-04 

O.OOiO.OO 

14.5±1.5 

HIP 74752 

1.63 

401±7 

138±5 

6.23e-04 

O.OOiO.OO 

4.1±0.45 

HIP 74959 

1.52 

68.6±0.6 

15.9±0.8 

2.38e-02 

l.OOiO.OO 

21.8±2.2 

HIP 75491 

1.52 

125±0.2 

5.12±0.05 

1.41e-03 

O.OOiO.OO 

13.5±1.3 

HIP 75509 

1.58 

251±1 

5.13±0.07 

4.37e-05 

O.OOiO.OO 

3.09±0.5 






Table 3—Continued 


HIP ID 

xl 

Tgr 

(K) 

ttgr 

(^m) 

mass 

(Afmoon) 

fo 

rgr 

(AU) 

HIP 76395 

1.39 

306±1 

9.53±0.07 

1.13e-04 

O.OOiO.OO 

5.32±0.73 

*HIP 77432 

2.17 

228±3 

8.75±0.4 

2.19e-05 

O.OOiO.OO 

1.62±0.63 

HIP 77520 

1.36 

243±4 

5.63±0.26 

5.68e-06 

O.OOiO.OO 

1.44±0.19 

HIP 77523 

1.02 

268±3 

8.96±0.22 

9.43e-05 

O.OOiO.OO 

6.56±1 

HIP 779II 

0.04 

112±1 

29.4±1 

1.02e-01 

l.OOiO.OO 

38.5±5.6 

HIP 78663 

0.95 

67±0.9 

2.23±0.41 

6.40e-03 

O.OOiO.OO 

91.8±9.1 

HIP 78756 

0.68 

204±1 

5.26±0.13 

5.96e-05 

O.OOiO.OO 

9.55±1 

*HIP 78977 

3.18 

1500±0 

2690±10 

8.00e-04 

O.OOiO.OO 

0.0781±0.01 

HIP 79400 

1.24 

238±2 

9.62±0.2 

7.78e-05 

O.OOiO.OO 

3.95±0.63 

HIP 79439 

0.69 

236±0.3 

56.2±1 

8.27e-04 

l.OOiO.OO 

15.5±2.2 

HIP 79710 

1.08 

245±1 

4.63±0.09 

3.30e-05 

O.OOiO.OO 

2.54±0.91 

*HIP 79742 

2.48 

93.6±0.3 

10.5±0.1 

2.33e-02 

l.OOiO.OO 

14±2.1 

HIP 79878 

1.40 

203±1 

9.53±0.22 

2.43e-04 

O.OOiO.OO 

15.9±3.6 

*HIP 79977 

2.54 

102±0.2 

ll.liO.l 

4.44e-02 

l.OOiO.OO 

11.5±2.7 

HIP 80024 

1.70 

169±0.1 

10±0.04 

1.81e-03 

O.OOiO.OO 

37.6±4.9 

HIP 82218 

1.70 

154±2 

9.54±0.47 

2.74e-04 

O.OOiO.OO 

5.3±0.71 


*Fits with formal > 2 but appear to be fit well with a single-grain model by visual 
inspection. 







Table 4. Fits to two-grain model 


HIP ID 

xl 

Ti 

ai 

massi 

fo,l 

n 

T2 

02 

mass2 

fo,2 

r2 



(K) 

(/zm) 

(ikfrnoon) 


(AU) 

(K) 

(^m) 

(Afrnoon) 


(AU) 

*HIP 53524 

2.76 

64.7±0.2 

20±0.3 

2.24e-01 

l.OOiO.OO 

42.2±0.2 

413±23 

3.32±0.56 

2.86e-06 

0.37±0.30 

1.24±0.14 

*HIP 55188 

2.20 

128±1 

13.1±0.3 

1.03e-02 

l.OOiO.OO 

13.2±0.7 

535±30 

3.91±0.62 

4.95e-06 

0.65±0.09 

1.06±0.13 

*HIP 58220 

2.20 

169±1 

6.61±0.02 

1.03e-04 

l.OOiO.OO 

4.1±0.39 

1500±0 

14.2±0.2 

1.72e-06 

l.OOiO.OO 

0.068±0.006 

HIP 58720 

1.23 

129±1 

52.2±1.8 

2.89e-02 

0.32±0.14 

35.5±4.1 

468±26 

6.87±1 

1.49e-05 

0.66±0.22 

2.94±0.47 

HIP 59502 

1.20 

104±0.4 

3.62±0.002 

2.67e-03 

O.OOiO.OO 

31.9±4.1 

368±5 

7.96±0.38 

4.21e-05 

0.95±0.14 

1.98±0.26 

*HIP 59898 

2.73 

153±2 

7.46±0.2 

1.48e-03 

0.53±0.09 

15.4±2.8 

524±51 

5.09±1 

7.52e-06 

0.03±0.18 

1.78±0.47 

tHIP 60183 

1.34 

150±3 

16.1±0.7 

2.65e-03 

0.35±0.12 

18±1.5 

314±56 

11±9.3 

2.21e-05 

l.OOiO.OO 

4.42±1 

HIP 60561 

0.78 

134±3 

25.8±1.6 

2.24e-03 

l.OOiO.OO 

16±2.8 

413±24 

19.2±3.8 

2.48e-05 

0.76±0.36 

1.86±0.38 

HIP 61684 

0.70 

106±1 

14.4±0.5 

6.85e-03 

l.OOiO.OO 

15.4±1.1 

267±2 

6.53±0.02 

1.72e-05 

O.OOiO.OO 

2.41±0.18 

*HIP 61782 

2.86 

115±0.3 

13.4±0.2 

2.93e-02 

0.73±0.08 

15.2±1.2 

1500±0 

4.1±0.12 

2.12e-07 

O.OOiO.OO 

0.205±0,02 

*HIP 62657 

6.13 

35.9±0.5 

0.539±0.05 

5.58e-h00 

O.OOiO.OO 

958±172 

119±1 

8.03±0.54 

2.08e-03 

0.47±0.22 

7.34±Do 

0.951±0^ 

HIP 63439 

1.14 

69.5±0.9 

15.6±1.4 

3.47e-02 

l.OOiO.OO 

24.3±2.9 

405±35 

2.69±0.73 

1.50e-06 

0.66±0.36 

HIP 63839 

1.14 

192±3 

5.67±0.16 

3.20e-04 

0.17±0.10 

7.82±0.98 

753±140 

4.92±2 

2.26e-06 

0.47±0.14 

0.817±0.32 

tHIP 63886 

0.85 

69.2±2 

15.7±2.9 

3.72e-02 

l.OOiO.OO 

28.9±7.9 

292±21 

6.23±1 

8.30e-06 

O.OOiO.OO 

1.66±0.5 

tHIP 64184 

1.07 

99±1.3 

15.6±0.9 

3.51e-02 

0.10±0.20 

11.3±0.5 

140±2 

2.93±0.26 

5.07e-04 

0.15±0.16 

7.38±0.3 

HIP 65875 

1.45 

107±1 

12.5±0.4 

3.30e-02 

l.OOiO.OO 

13.2±1.4 

208±7 

1.86±0.39 

4.43e-05 

O.OOiO.OO 

4.59±0.55 

HIP 65965 

1.20 

108±1 

16.5±0.7 

1.39e-02 

0.90±0.15 

27.4±3.2 

367±12 

17.1±2 

5.50e-05 

0.66±0.32 

2.76±0.36 

HIP 66068 

1.51 

146±4 

23.8±1.9 

4.08e-03 

0.65±0.41 

15.9±4.8 

378±20 

3.91±0.81 

1.34e-05 

0.67±0.21 

2.65±0.84 

HIP 66566 

0.98 

153±2 

14.6±0.6 

1.73e-03 

0.52±0.20 

10.8±1 

428±41 

4.08±1 

4.22e-06 

O.OOiO.OO 

1.67±0.35 

tpiP 67230 

0.27 

161±3 

5.06±0.4 

4.03e-04 

l.OOiO.OO 

7.55±0.84 

1500±0 

2.13±0.2 

7.68e-08 

O.OOiO.OO 

0.213±0.02 

HIP 67970 

0.77 

139±2 

7.44±0.31 

6.67e-04 

0.65±0.30 

6.4±0.69 

1020±240 

6.31±2 

4.25e-07 

0.82±0.24 

0.188±0.09 

tHIP 68080 

1.02 

157±5 

23.9±2.4 

2.90e-03 

0.49±0.35 

22.4±3.4 

682±149 

17.8±9.7 

9.61e-06 

l.OOiO.OO 

1.47±0.68 

tHIP 70149 

0.85 

108±0.5 

7.64±0.01 

7.01e-04 

O.OOiO.OO 

10.1±1.4 

268±8 

4±0.54 

6.75e-06 

O.OOiO.OO 

1.66±0.26 

HIP 70441 

1.87 

69.7±0.8 

20.9±1.4 

5.23e-02 

l.OOiO.OO 

46.3±5.8 

268±5 

3.9±0.26 

1.50e-05 

O.OOiO.OO 

3.27±0.42 

HIP 71271 

1.13 

84.6±1.4 

21.9±1.9 

3.69e-02 

l.OOiO.OO 

47±5.9 

311±14 

3.96±0.66 

1.07e-05 

O.OOiO.OO 

3.74±0.56 

*HIP 73145 

2.24 

80.4±0.3 

39.1±0.7 

5.60e-01 

l.OOiO.OO 

35.1±7.5 

214±1 

4.5±0.16 

2.77e-04 

O.OOiO.OO 

4.65±1 

HIP 73990 

1.63 

166±2 

9.91±0.39 

8.26e-04 

1.00±0.19 

6.62±0.86 

1440±890 

2.43±2 

1.42e-07 

0.46±0.23 

0.173±0.22 






Table 4—Continued 


HIP ID 

xl 

Ti 

(K) 

ai 

(/xm) 

massi 

(Afmoon) 

/o,l 

ri 

(AU) 

T2 

(K) 

02 

(^m) 

mass2 

(Afinoon) 

fo,2 

r2 

(AU) 

ttfflP 74499 

0.75 

83.7±0.8 

12±0.7 

7.45e-03 

l.OOiO.OO 

13.1±1.6 

286±3 

3.98±0.01 

3.21e-06 

O.OOiO.OO 

1.14±0.14 

HIP 75077 

0.39 

63.5±0.9 

30.4±2.6 

4.48e-02 

O.OOiO.OO 

79.4±14 

291±14 

4.08±0.63 

5.61e-06 

O.OOiO.OO 

3.92±0.8 

^HIP 75210 

0.81 

151±4 

139±10 

2.51e-02 

0.42±0.34 

26.4±3.2 

477±59 

4.91±1 

6.10e-06 

l.OOiO.OO 

2.73±0.74 

*HIP 76310 

2.19 

116±0.4 

17±0.2 

5.46e-02 

l.OOiO.OO 

52.9±5.2 

361±12 

10.8±1.3 

6.72e-05 

0.40±0.27 

6.09±0.71 

Imp 77315 

0.52 

153±4 

22.9±1.5 

5.16e-03 

0.26±0.28 

18.4±1.8 

440±48 

7.23±2 

1.63e-05 

0.76±0.29 

2.54±0.59 

HIP 77317 

0.88 

97.3±1 

54.2±3.9 

6.14e-02 

0.96±0.32 

35±7.8 

486±35 

12.8±2.6 

6.66e-06 

0.42±0.51 

1.56±0.41 

HIP 78043 

0.75 

78.2±1.7 

20.7±2.5 

2.88e-02 

l.OOiO.OO 

23.1±4 

307±13 

6.05±0.97 

9.96e-06 

O.OOiO.OO 

1.47±0.27 

HIP 78641 

1.54 

131±1 

10.5±0.4 

2.47e-03 

0.83±0.13 

9.12±1 

366±17 

2.32±0.41 

4.19e-06 

0.06±0.12 

1.5±0.32 

HIP 78996 

1.72 

175±1 

9.93±0.01 

3.61e-04 

l.OOiO.OO 

14.1±2 

1500±0 

4.96±0.08 

5.40e-07 

0.71±0.09 

0.306±0.04 

HIP 79156 

1.45 

195±3 

9.06±0.39 

2.55e-04 

l.OOiO.OO 

18.9±2.9 

1500±0 

5.29±0.11 

1.19e-06 

0.89±0.09 

0.476±0.07 

HIP 79410 

1.20 

114±4 

8.02±0.91 

1.37e-03 

O.OOiO.OO 

57.5±9.9 

387±9 

6.71±0.51 

3.65e-05 

l.OOiO.OO 

5.62±0.a§ 

HIP 79516 

1.77 

72.1±1.5 

24±2.6 

2.51e-01 

l.OOiO.OO 

24.9±4.5 

128±3 

3.32±0.59 

5.67e-04 

0.11±0.15 

9.62±1 

*HIP 79631 

2.33 

76.8±0.8 

35±0.2 

2.90e-01 

0.09±0.36 

66±10 

205±4 

7.4±0.78 

3.51e-04 

O.OOiO.OO 

8.78±1 

t*HIP 80088 

2.86 

51.7±0.4 

1.55±0.11 

2.45e-01 

l.OOiO.OO 

1000±140 

245±3 

3.4±0.17 

3.10e-05 

O.OOiO.OO 

8.37±1 

Imp 80320 

0.84 

50.6±0.9 

0.979±0.16 

8.23e-02 

l.OOiO.OO 

359±68 

337±18 

2.34±0.43 

2.74e-06 

0.23±0.23 

1.21±0.26 

*HIP 80897 

2.90 

108±1 

93.5±2.2 

1.75e-01 

l.OOiO.OO 

30.2±5.1 

481±16 

36.3±3.5 

6.89e-05 

l.OOiO.OO 

1.54±0.28 

HIP 82154 

0.86 

98.3±2.5 

22.3±2.2 

4.74e-02 

O.OOiO.OO 

75.7±19 

257±8 

6.35±0.88 

1.93e-04 

l.OOiO.OO 

11.1±2.9 


*Fits with formal > 2 but appear to be fit well with a two-grain model by visual inspection. 

^Formally well-ht by a single-grain model, but but appear to be better fit by a two-grain model by visual inspection. 
TRS spectra well-fit by a two-grain model, but does not fit the 70 /xm MIPS photometric value. 







Table 5. Objects not fit well by single or two-grain models 


HIP ID 

xl 

Ti 

Oi 

massi 

fo,l 

n 

T2 

02 

mass2 

fo,2 

r2 



(K) 

(/an) 

(Aifrnoon) 


(AU) 

(K) 

(^m) 

(Afixioon) 


(AU) 

HIP 63975 












one-grain 

78.77 

478±0.2 

4.54±0.004 

8.45e-04 

I.OOiO.OO 

2.51±0.31 






two-grain 
HIP 79288 

49.00 

101±0.06 

2.62±0.0003 

8.06e-02 

0.49±0.01 

74.6±9.1 

lllOilO 

2.5±0.02 

6.68e-05 

l.OOiO.OO 

0.65±0.08 

one-grain 

793.62 

361±0.1 

4.38±0.01 

5.49e-04 

O.OOiO.OO 

1.35±0.2 






two-grain 

138.68 

lOOiO.l 

2.04±0.02 

4.7Ie-02 

0.51±0.01 

33.9±5.1 

822±4 

3.56±0.03 

4.56e-05 

0.48±0.00 

0.418±0.06 
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Fig. 5.— Best fitting single grain and two-grain models to HIP 56673, with the photosphere 
£t to the 5-6 micron region of the IRS spectrum. We find that a single component model 
best fits this spectrum. 


HIP 56673* 

1- grain (0.9) 

2- grain (0.9) 



This count does not include HIP 56673 or HIP 78977. As discussed previously, these 
sources have excesses similar to Rayleigh-Jeans profiles, consistent with either a photosphere 
mismatch or a hot dust component. If the photosphere is scaled to match the 5-6 micron 
region of the IRS spectrum, then HIP 78977 is a non-excess source, and HIP 56673 has a 
marginal excess that can be modeled well as a single belt. In Figure [5l we show the best fit 
to HIP 56673, after fitting the photosphere to the IRS spectrum instead of calibrating the 
photosphere to optical and near-IR photometry. 


Two of our targets, HIP 63975 (HD 113766) and HIP 79288 (HD 145263) have high 
amounts of excess and show evidence of silicate features that cannot be adequately repro¬ 
duced with the simple models used in this paper. The mineralogy of HIP 63975 was analyzed 
in detail in iLisse et al.l (120081) . and was found to consist of both amorphous and crystalline 
silicates, as well as Fe-rich sulfides, amorphous carbon, and water ice. HIP 79288 has a 
similarly complex composition, and requires silica in addition to other species to adequately 
model it (Lisse, et ah, in prep). 


HIP 74499 and HIP 80088 both appear to require at at least two components to fit 
their IRS spectra: a large cold component to fit the 20-40 micron region, and a small hot 
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component to fit the 10 micron region. However, in both these cases, the best £t to the IRS 
spectrum does not produce enough emission at 70 micron to match the MIPS photometry. It 
is likely that these two systems have an additional third cold component that would account 
for the missing 70 micron excess. Both these sources have been classihed as two-belt systems 
for further analysis. 


Recently, a debris belt has been imaged around one of our targets, HIP 64995 (HD 
115600), with a semi-major axis of ~ 48 AU. (jCurrie et al.ll2015l) . Our modeling predicts 
a single belt of debris at a temperature of 114 K, or ~ 11 AU. This belt woul d be interio r 


to the inner working angle of the coronagraphic image of the belt detected by ICurrie et ah 
(120151) . If the outer belt is an analog of the Kuiper belt, then the belt predicted by the mid- 
infrared spectroscopy would be an asteroid belt analog. This example illustrates how mid-IR 
spectral analysis is complementary to coronagraphic imaging for studying the structure of 
debris disks. 


In Figure [6l we show the masses of the dust belts in the single-belt (black) and two-belt 
(red/blue) debris disk systems. For the two-belt systems, the hot and cold belts are colored 
red and blue, respectively. Dust mass appears to be inversely correlated to temperature. 
This is likely a selection effect, because more cold material must exist in order for it to 
contribute signihcantly to the infrared excess. That is, a lower temperture blackbody emits 
less radiation total than a higher temperature one, holding the emitting surface area constant. 
There does not appear to be a trend in measured dust mass versus stellar mass. This implies 
that the trend that Lir/L* decreases with increasing stellar mass (see Figure E]) is best 
explained by differences in dust temperature (i.e. distance from the system primary) rather 
than differences in measured dust mass. 


In Figure [TJ we summarize the results of our models in terms of the temperatures 
and radial distance of the grain populations. In general, the temperatures of the single-belt 
models are intermediate between the hot and cold components of the two-belt models. There 
appears to be more scatter in the hot components for lower stellar masses, which is more 
easily seen in the cumulative distributions in temperature after separating them out into low 
mass (< 1.5 Mq, dashed lines) and high mass (> 1.6 Mq, solid lines) stars. 

The same discrepancy between high-mass and low-mass stars appears in the distribution 
of radii of the belts, as shown in the right panel of Figure [71 The belt radii are an average 
distance for the olivine and pyroxene components of the grains. The model spectra assume 
that the olivine and pyroxene components of the dust all have the same equilibrium temper¬ 
atures. The grain properties are similar enough between olivine and pyroxene that grains 
of the same temperature are effectively co-located as well, and the radius derived for pure 
olivine is nearly the same as that of pure pyroxene. Since in many cases the composition is 
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Fig. 6.— Dust belt masses versus dust temperature (top) and stellar mass (bottom), for the 
single grain population disks (black) and the two-grain population disks (red/blue). In the 
two-grain population model, the cold component is shown in blue while the hot component 
is shown in red. Objects belonging to LCC, UCL, and USco are labeled by asterisks, squares, 
and diamonds, respectively. 
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Fig. 7.— Temperature and radial distribution of grain populations in single belt (black) 
and two-belt (blue/red) systems. For the two-belt systems, the hot component is plotted 
in red, the cold component in blue. Top row: Tgr (left) and r (right) versus stellar mass. 
Objects belonging to LCC, UCL, and USco are labeled by asterisks, squares, and diamonds, 
respectively. Bottom row: Cumulative temperature distribution (left) and radial distribution 
(right) for the dust components of best £t disk models. The solid lines show stars with mass 
> 1.6 Mq, while the dashed lines show stars with mass < 1.5 Mq. 













43 


either pure pyroxene or pure olivine, the difference in the distances of the two components 
is insignihcant. 


The grain temperatures of dust around high mass stars are sys t emati cally higher than 
around low mass stars. This agrees with the finding in IChen et al.l fl2014j) . In Table [6l we 
show the results of the Kolmogorov-Smirnov (KS) and Anderson-Darling (AD) statistics for 
testing the similarity of the grain temperatures between high and low mass stars. We consider 
the single belt systems, the cold component of the two-belt syste ms, and the hot component 
of the two-belt systems separately. The Anderson-Darling test flAnderson fc Darlind 1195411 
is arguably a more sensitive statistic than the KS test because it gives more weight to the 
the tails of the distributions. To test whether or not two samples Xi,, Xn and Yi,..., 
are drawn from the same populations, we can use the two-sample Anderson-Darling statistic, 
given by 


= 


mn 


N-l 

Y 

2=1 


{MiN - nif 
i{N -i) 


(17) 


where N = m + n and M j is the numb er of X’s less than or equal to the Ah smallest element 
in the combined sample fjPettittlll976l) . 


The KS test for the single belt systems and the cold component of the two belt systems 
give a 10% and 4% or less probability, respectively, of being drawn from the same popula¬ 
tion, while the probability of the hot components being similar is 37%. The AD test gives 
comparable results: < 5% probabilities for the single belt systems and cold components, but 
higher probability of similarity for the hot component (> 15%). A potential explanation for 
this systematic shift is that higher mass stars heat the dust more than lower mass stars do, 
but that is not is fully borne out when we examine the distributions with respect to dust 
distance. 


It appears likely that the cold grain populations for both the high mass and low mass 
stars are at similar distances: the KS test gives a probability of 65% that they are drawn 
from the same population. However, both the single belt systems and the hot component of 
the two-belt systems are systematically offset to larger distances in the high mass systems, 
despite being hotter. The single belt systems have a 22% probability of being similar, while 
the probability for the hot component of the two belt systems being similar is 12%. Again, 
the AD test gives similar results: it is unlikely that the single belts and hot components 
of low mass stars are drawn from the same population as the high mass stars (< 10%). 
However, it is much more likely that the cold components are similar (> 15%). 

The trend in Ljr/L^ seen in Figure [3] also supports the idea that lower mass stars have 
dust at closer radii. This is because the infrared excess luminosity, which is attributed to 
the dust, scales as L/r oc where and are the total mass and temperature of 
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the dust, respectively. The temperature of the dust is determined by stellar illumination, as 
oc L^/d? where d is the stellocentric distance. This gives Ljr/L^ oc Md/R^. If Md stays 
constant, then the trend that lower mass stars have higher Ljji/L^ implies that the disks in 
lower mass stars are more compact. 


This implies that that the low mass stars retain close-in dust more readily than high 
mass stars, suggesting that debris disks in high mass stars evolve faster than low mass stars, 
and that this evolution occurs inside out. One explanation for thi s is that debris disks evolv e 
faster in high mass stars because the dynamical times are shorter (IKenvon fc Bromlevll2008l) . 
Another possibility is that the higher mass stars (F-type and earlier) evolve onto the main 
sequence sooner. By 15-17 Myr, the ages of LCC and UCL, these stars are already on the 
main sequence. The ignition of hydrogen burning in these stars could enable the clearing 
out of inner dust belts at < 1 AU. The retention of dust by low-mass stars could also simply 


be the result of more efficient grain blow-out by more massive stars, as indicated by their 
larger values of amin [Eq. ([7])]. Still another possibility is that the initial protoplanetary gas 
disk differs between high mass and low mass stars. 


5. Discussion 


The high degree of scatter in the distances of the dust belts indicates that they are 
not likely connected to intrinsic properties of the primordial disks from which they arose. 
Primordial disks generally have continuous radial distributions of material rather than belts. 
Material could potentially pile up at pressure maxima, for example, but the origin of the 
pressure maxima depends on stellar properties such as effective temperature and luminosity, 
which should be relatively stable for a given stellar mass. An example of this phenomenon 
i s the T Tauri object HL Tau, w hich has been seen to have several gaps in ALMA imagery 
flALMA Partnership et al.N2015l) . Although it is possible that these gaps were created by 
plane ts, locations of th e gaps also appear to be co-incident with condensation fronts in the 
disk fjZhang et al.ll2015l) . The locations of the fronts are based on the predicted temperature 
profile of the disk, which in turn dep ends on the heating of the disk from stellar irradiation 
(see e.g. 


Kenvon fc HaTtma.nnlll98 


r 


Alternatively, dust locations could depend on the formation of planets, which is highly 
stochastic in regard to stellocentric distances. The varying diameters of the inner holes seen 
in transitional disks are often attributed to planet formation for this reason. Planets could 
also explain the origin of the dust belts after the dissipation of the primordial gas, since 
planets could be responsible for shepherding the parent bodies that produce the dust. 
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Stellar companions could also affect the dust belts in our disks. A high fraction of stars 
in our sample have identified stellar companions. A few stars have directly imaged substellar 
companions that may be distant planets. Since planets and binary companions are likely to 
play key roles in the sculpting of debris disks, it is important to put the properties of dust 
belts into context with the presence of binarity and the presence of planets. 


5.1. Binary Stars 


A number of s tars in our sample have been identihed to have binary companions 


(Janson et ah 

2013; 

Chen et ah 

2012: 


Kouwenhoven et al.ll2005L l2007|l . The projected dis¬ 


tances of these binary companions are tabulated in Table [71 The distances are calculated 
from the angular separati ons from the above references and using Hipparcos stellar distances 
from Ivan LeeuwenI (1200711 . 


Stellar companions will truncate a circumstellar disk through tidal interactions. For a 
debris disk, where gas has little to no dynamical effect, the disk truncation radius can be 
estimated from the last stable orbit. Assuming a circular orbit, the outermost radius of a 
circumstellar disk allowed by a binary companion can be expressed as 


flint = (0.464 — 0.380/i)a 
and the inner edge of a circumbinary disk is 


(18) 


flext = (1-60 -|- 4.12p,)a 


(19) 


flHolman fc WiegertlllOOOll . where fi is the ratio of the mass of the binary companion to the 
total masses of the two stars, and a is the semi-major axis of the binary orbit. 


In Figure [HI we show the outer location of the best-fit dust belts versus the binary 
separation for those objects with binary companions. If the dust is best fit with a single 
grain population, then that distance is used. For the remaining systems, the location of 
the outer belt in the two-grain fit is used. We also indicate the disk truncation radius for 
binaries with mass ratio 0.5 (equal mass binary) and 0.1 assuming a circular orbit. 


In all cases, the dust is located interior to the binary separation, so the binary companion 
must have truncated the disks in all these systems. A few dust belts appear close to the 
truncation radius of a /x = 0.1 binary companion, however the projected binary separation is 
affected by the inclination of the system and eccentricity of the orbit, so the true truncation 
radius might be higher. No circumbinary dust belts in our sample of stars has been identified. 
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Table 6: Kolmogorov-Smirnov and Anderson-Darling tests for grain populations 

population KS test AD test 

D p(D) Al„ p(Ay 


Temperatures 


... single belt 

0.33 

0.097 

2.5 

< 0.05 

... cold component 

0.45 

0.045 

2.6 

< 0.05 

... hot component 
Distances 

0.30 

0.37 

1.5 

> 0.15 

... single belt 

0.28 

0.22 

2.1 

< 0.10 

... cold component 

0.24 

0.65 

0.93 

> 0.15 


... hot component 0.39 0.12 

^Based on tabulations of probabilities in 


1.9 < 0.10 


StephensI (^19741) . 


Table 7. Projected Distances of Binary Companions 


HIP 

HD name 

ai (AU) 

02 (AU) 

03 (AU) 

53524 

HD 95086 

440.3 



58220 

HD 103703 

75.17 



58528 

HD 104231 

492.8 



59502 

HD 106036 

296.1 

911.4 


59693 

HD 106389 

58.49 



63236 

HD 112383 

837.2 

1326. 


63839 

HD 113457 

427.4 

595.4 

627.2 

65965 

HD 117484 

1514. 



68080 

HD 121336 

268.5 

3902. 


69291 

HD 123889 

195.4 



77315 

HD 140817 

100.1 

5552. 


77317 

HD 140840 

4126. 



77520 

HD 141254 

223.5 



78756 

HD 143939 

1249. 



79400 

HD 145357 

734.2 



79631 

HD 145880 

376.0 

1132. 


82154 

HD 151109 

1869. 
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Fig. 8.— Outer dust belt location versus projected binary separation. The solid line marks 
the 1:1 line, so circumstellar disks lie below the solid line while circumbinary disks lie above 
the solid line. The dotted and dashed lines mark the disk truncation radii for an equal-mass 
binary and a /r = 0.1 binary, respectively. Crosses indicate single belt systems, and asterisks 
indicate the outer belt of two-belt systems. Crosses indicate systems not well-modeled by 
either single or two-belt systems, but using the best htting two-belt model regardless. 


5.2. Stars with planets 


Three of our sources have bee n identihed to be planet hosts. HIP 53524 (HD 95086) 
is a binary star fjChen et al.l l2012l ) around which a ~ 5 Mt„p planet has been detected at 
a projected separation of 56 AU by direct imaging (iRameau et al.l 120131) . HIP 59960 (HD 
106906) hosts a very distant planetar y-mass companion at a projected separation of ~ 650 
AU, also detected by direct imaging f Bailev et ah 2014). HIP 73990 (HD 1 33803) has at 


least two planets in its system, at projected separations of 20 and 32 AU flHinklev et ah 
2015h . 


HIP 53524 is best £t by a two-belt model, though its relatively high reduced chi-square 
value {xl = 2.8) suggests that a few more addition parameters, such as hnite belt widths or 
the addition of crystalline silicates, could improve the model. It hosts a stellar companion 
at 440 AU in addition to the imaged planet. The positions of the belts are 42 AU (65 K) 
and 1.2 AU (413 K), both interior to the projected position of the planet. Since the outer 
belt is close to the projected planet distance of 56 AU, it is likely that it is sculpted by the 
planet. 


Herschel images of HIP 53524 marginally resolve the disk, suggesting that the systems 
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is surrounded by a halo that extends to around 800 AU fjSu et ahl 1201411 . That analysis 
included Spitzer/MIPS data and proposed that in addition to the halo, the systems consists 
of a warm belt at 175 K and a cold belt at 55 K. The wavelength coverage of Spitzer does 
not provide much sensitivity to emission from the cold halo, but the 55 K belt is consistent 
with the model for HIP 53524 presented here. Our model predicts a hotter inner belt, which 
may be a result of not including the third coldest dust distribution that is inferred from the 
Herschel data. 

In our analysis, HIP 59960 is well-h t by a single dust belt at 12 AU, and its reported 
planetary companion flBailev et al.ll2014j) is so distant that it does not interact with the dust. 
It is possible that additional unseen companions exist between the belt and the planetary 
companion and that those companions could shephered the dust. 

HIP 73990 is best £t to a two-belt model with grain temperatures of 1440 K and 166 
K, corresponding to distances of 0.17 and 6.6 AU, respectively. Assuming that planets are 
able to cle ar material out to t he 2:1 mean motion resonance, then the inner planet imaged 
at 20 AU f Hinklev et ah 2015) should truncate the debris disk to 12.6 AU. Both dust belts 
detected from our spectroscopic modeling are interior to this distance. 

Additional objects in ScoCen that have detected sub-stellar companions include HIP 
78530, which exhibits no i nfrared excess; GS C 06214-00210, an accreting T Tauri star; and 
IRXS J160929.1-210524. f Bailev et ah 2013 1. These objects were not included in our study 
because they are pre-main sequence stars and their disks are protoplanetary in nature. 

Although imaging planets around debris disk systems can help us understand the role 
that planets play in sculpting debris disks, few such planet images exist. On the other hand, 
many debris disks are well-studied. We can turn the question around, then, and ask what can 
be learned about planet formation from debris disks. Dynamical interactions with planets 
should sculpt and shephered debris disks. Therefore, debris disks with structure, such as 
gaps, can imply the presence of unseen planets. 

Numerical simulations have shown that a companion orbiting in a disk can create gaps 


via p lanetesimal scattering in overlapping resonances (e.g., Roques et al.lll994 : Lecavelier des Etangs et al 


1996ll. The width of the gap is r e lated to the rnass of the companion by a power law (e.g.. 


QuillenI 120061: IChiang et al. 


20091: iRodigas et ahl 1201411 . t he parameters of which d epend on 


the age of the system and the optical depth of the disk flNesvold fc Kuchnen 1201411 . 


Following the procedure described in iNesvold fc Knchnerl (1201411 . we analyzed the two- 
belt systems in our sample to place upper limits on the mass of a possible single perturbing 
companion in each system. In each case, we assumed that a single body on a circular orbit 
equidistant in log semimajor axis between the two dust bands has cleared the gap between 
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the bands. We used as a proxy for the face-on optical depth of each disk. 


Nesvold &: Kuchnerl (120141) found that the largest gap size that a single body on a 


circular orbit can create has a full width of Ar/r = 1.6. Larger bodies tend to stir the disk; 
destroying it and widening it, while roughly preserving the gap edge near the location of the 
2:1 mean motion resonance. Nine disks in our sample had gaps narrower than this maximum 
width. Table ([8]) summarizes inferred companion masses and semimajor axes for these nine 
disks. 


Five of these systems, HIP 61684, HIP 66068, HIP 78641, HIP 65875, and HIP 79516, 
have small enough gaps that they require a single perturbing body whose mass is in the range 
of planet masses. The other disks require either planets on eccentric orbits, companions in 
the brown-dwarf mass range or multiple planets. The inferred companions in these systems 
are all located within 0.2 arcseconds of their host stars, probably too close to detect directly 
with today’s instruments. However, future observatories, such as WFIRST-AFTA, may be 
able to resolve sub-stellar companions at angular separations of less than 0.2". In addition, 
the next generation of large aperture ground-based telescopes, such as GMT, TMT, or E- 
ELT, could have the necessary resolving power and inner working angle. ATLAST, currently 
a NASA strategic mission concept study, could also detect these planets. 

HIP 82154 also has a binary companion at a projected distance of 1869 AU (see Table 
IH]). This companion is probably too distant to have any dynamical effect on the dust belts 
or any planet located between them. 

Six out of the nine inferred companions listed in Table (jS]) have estimated semimajor 
axes within 5 ± 3 AU, close to Jupiter’s semi-major axis of 5.2 AU. This lends credence to 
the idea that giant planets preferentially form in the 5-10 AU range. This is also consistent 


Table 8: Masses and semimajor axes of companions inferred from the disk gaps. 

HIP Companion Semimajor Axis (AU) Max Companion Mass {Mjup) 


60561 

5.4 

28.8 

61684 

6.4 

12.2 

65089 

22.9 

33.1 

66068 

5.9 

10.4 

66566 

4.2 

15.3 

78641 

3.7 

9.8 

82154 

28.7 

42.9 

65875 

7.7 

1.1 

79516 

15.3 

0.8 
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1 10 

Semimajor Axis (AU) 


Fig. 9.— Maximum companion mass vs. estimated companion semimajor axis for the nine 
companions inferred from the disk gaps. The bars indicate the widths of the gaps between 
belts. 


with statistics of exoplanets showing that, apart from hot Jupiters, the number of exoplanets 
increases toward orbital separations larger than 1 AU. 
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6. Conclusions 

Constraining our study of debris disks to those in ScoCen gives us the advantage of 
examining a cohort of debris disks of similarly young age. Within this sample, we hnd 
evidence of mass-dependent evolution of the hot dust. In particular, in systems with two 
belts, low mass stars have closer inner belts than high mass stars. This implies that high 
mass stars have less hot dust than low mass stars. This could be related to the faster 
evolution times of the higher mass stars, which results in higher mass stars reaching the 
main sequence sooner than the less massive stars. Then the lack of hot dust could be 
explained by the clearing of the dust from the inside out as the star evolves. 

We explored how stellar and sub-stellar companions could sculpt debris disks. Many of 
the objects in our sample have known binary companions. We hnd that the dust distances 
from our models predict circumstellar disks rather than circumbinary disks. However, be¬ 
cause the binaries are generally wide, and our data is limited to infrared wavelengths, our 
observations are not sensitive to any thermal emission that might come from circumbinary 
disks. The dust distances are consistent with disk truncation at outer radii by the binary 
companion. 

Two of our objects host known planets. These planets have been detected by direct 
imaging, so they are distant planets outside the dust disks. The fact that the planets that 
have been discovered in ScoCen also host debris disks suggest that debris disks and planet 
formation are correlated. Planets can also sculpt debris disks, and we use this fact to consider 
the possibility that our disks might host planets. 

The locations of dust belts can put constraints on the locations of planets since they 
should clear out gaps in the disks. The two-belt systems found in our study could be the 
result of one or several planets carving out gaps. If the distance ratios are small, then the 
location and mass of the potential planet can be narrowly constrained. These systems are 
particularly good targets for follow-up planet searches. 

Wider gaps could be created by multiple planets or eccentric orbits. These are also good 
targets for follow-up, although predictions about planet propeties are less well-constrained. 


This work is based on observations made with the Spitzer Space Telescope, which is 
operated by JPL/Caltech under a contract with NASA. Support for this work was provided 
by NASA through an award issued by JPL/Caltech. HJ-C acknowledges support from NASA 
grant NNX12AD43G. 
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